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PREFACE 


This  Report  is  one  in  a  series  entitled  '  Effects  of  Strong  Explosions" 
and  intended  to  provide  an  overview  of  current  Soviet  research  activities  in 
that  field.  It  falls  within  the  scope  of  a  continuing  program,  sponsored  by 
the  Defense  Advanced  Research  Projects  Agency,  which  undertakes  the 

systematic  coverage  of  selected  areas  of  Soviet  scientific  and  technological 
literature . 


Thp  effects  of  strong  explosions,  particularly  nuclear  explosions,  and 
their  simulation  under  laboratory  and  field  conditions  involve  a  very  broad 
range  of  phenomena  and  parameters  such  as  high  temperatures,  pressures, 
electromagnetic  energy  densities,  interaction  of  energy  with  materials, 
shock  waves  in  gases,  liquid,  and  solids,  and,  in  general,  events  occurring 
under  conditions  of  a  high  degree  of  mechanical,  thermal,  and  radiation 
stress.  The  series  of  Rand  Reports  of  which  this  is  the  fifth  to  be  published 
will  cover  a  variety  of  the  pertinent  topics  as  they  are  reflected  in  the  cur¬ 
rently  available  Soviet  literature,  by  supplying  abstracts  from  and  summaries 
of  the  scientific  and  technological  publications  of  the  USSR. 


The  material  in  the  present  Report  is  derived  exclusively  from  Soviet 
technical  writings  that  appeared  in  1971  and  the  early  part  of  1972,  The 
abstracts  here  included  were  prepared  by  Informatics,  Inc.  They  have  been 
arranged  by  the  compiler  according  to  fairly  general  subject  areas,  and  vary 

in  length  depending  on  the  pertinence  and  significance  of  the  individual 
papers . 

The  information  assembled  in  this  Report  covers  our  most  recent 
knowledge  of  Soviet  research  on  the  effects  of  explosions,  and  should 

prove  useful  to  Americ  in  scientists  working  in  this  and  related  fields  of 
research. 
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I.  SHOCK  WAVES  IN  GAS 


Aleksandrov,  V.  V.  Phase  plane  method  for 
solution  of  one-dimensional  problems  in 
radiative  gas  dynamics.  MZhiG,  no.  1,  1972, 

144-155. 

The  phase  plane  method  is  used  to  solve  the  problem  of  one-di¬ 
mensional  equilibrium  flow  of  an  inviscid,  radiating,  absorbing,  and  scattering 
gas.  For  gas  propagation  with  strong  radiation  interference,  the  usual 
numerical  solution  to  this  problem,  based  on  Peierls  equation  of  radiant 
emission  W as  a  function  of  optical  depth  t  ,  becomes  complicated  because  of 
boundary  layer  formation  behind  a  shock  wave.  For  that  reason  the  problem 
is  formulated  in  different  material  coordinates;  namely,  dimensionless  gas 
velocity  V  and  W,  in  the  (V,W)  or  "phase  plane".  In  the  presence  of  scatter¬ 
ing  and  strong  radiation  interference  with  gas  propagation,  such  an  approach 
is  preferable  to  using  the  Peierls  equation,  because  the  V  of  one-dimensional 
gas  flow  is  a  measure  of  both  kinetic  and  internal  gas  energies.  The  phase 
plane  thus  represents  an  energetic  space  with  the  energy  characteristics  of 
the  substance  and  radiation  as  coordinates.  The  function  W(V)  for  gray  body 
radiation  is  determined  by  the  nonlinear  integral  equation 

B  * 

•« 

,  (1) 

XEi{B\  j  A’[Tj,u>(T])]dr|  |  J  dl 


where  B  =  Bo/4  is  the  Boltzmann  constant,  0  is  the  dimensionless  gas 
temperature,  K  is  the  absorption  coefficient,  £ ,  Ej,  and  t]  are  values  from 
Peierls  equation  in  neutron  transport  theory.  Solution  of  (1)  gives  W(V). 

The  phase  plane  (V,W)  for  selective  radiation  has  an  infinite  dimension,  and 
the  problem  for  a  nonscattering  gas  is  therefore  formulated  by  the  equation 
of  energy 


X 


I 


|  *i|} 


(2) 


where  q  is  the  volumetric  rate  of  gas  energy  increase  due  to  radiation, 

Bs  =  2L _  ®°*  Kp,  K,  h,  y,  and  P  are  the  dimensionless  Planck  absorption 
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coefficient,  optical  constant,  enthalpy,  frequency,  and  Planck  function 
respectively.  ' 


Application  of  formulas  (1)  and  (2)  is  illustrated  by  examples 
of  a  strong  shock  wave  propagating  in  a  cold,  transparent,  nonscattering 
gas  and  a  shock  wave  propagating  in  a  radiating,  absorbing,  and  scattering 
gas  with  strong  radiation  interference.  In  the  first  example,  a  formula 
derived  from  (1)  for  the  0  of  a  gray  body  shows  that  the  gas  is  strongly 
cooled  by  radiation.  The  discontinuity  of  shock  wave  velocity  is  formulated 
m  the  second  example  for  a  perfect  gas  in  a  diffusion  approximation  of  the 

^  tHe  PhaSe  Plane*  is  shown  that  discontinuity 
exists  at  7  2  and  Mj>  1.5  or  >  2.  05  in  the  incident  flow. 


Simonov,  I.  V.  Diffraction  of  a  strong  shock 
wave  on  a  weakly-defined  wedge.  ZhPMTF, 
no.  6,  1971,  107-114. 


A  theoretical  study  is  presented  on  diffraction  of  a  strong 
wo  -  dimensional  steady  shock  wave  on  a  wedge  whose  properties  (initial 
density,  behavior  under  impact  load)  differ  only  slightly  from  those 

°l  medium  in  which  the  sho<*  wave  propagates.  A  linear  approximation 
of  diffraction  -  initiated  small  perturbations,  disregarding  differences  in 
strength  and  sound  velocity  c  in  the  medium  and  wedge,  is  justified  under 
conditions  of  high  pressure  (  ~  10^  atm)  and  temperature  behind  the  strong 
shock  wavefront.  The  validity  of  the  approximation  is  supported  by  the  fact 
t  at  the  effect  of  strength  characteristics  diminishes  with  increased  shock 
wave  intensity.  In  this  approximation,  the  unknown  perturbation  velocity 
components  u'  and  w'  and  perturbation  pressure  p'  satisfy  the  ordinary 
knear  equations  for  the  two-dimensional  flow  of  an  ideal  compressible  fluid. 
The  equations  are  given  in  the  form 


DP  — 


0u_ 

bx 


,  Du  =  °P,  Dw  =  °P  { 

Ox  '  0y  \ 


(1) 


where  u,  w,  and  p  are  dimensionless  variables.  Conditions  in  the  pertur¬ 
bation  region  of  the  shock  wavefront  are  described  by 


P  A  (f  —  y]')t  u  _  ftp*  w  _  at  x  =  k 

in  the  medium,  and 

p  =  A  (/  —  i /}')  +  Alt  u  =  IJp  +  Bx,  w  =  —Mf  at  x  —  k 
A  —  2x,/  A  _  ;.vr(2x— !)»'  —  y.ra  n  l+i 

1~k]I  -  5==-ir 


(2) 


(3) 


in  the  wedge.  Vo  and  Uo  in  (3)  are  specific  volume  and  mass  velocity, 
respectively,  behind  the  shock  wavefront,  and  V  is  the  specific  volume  of 
the  medium  ahead  of  the  front.  The  .problem  was  reduced  to  determining  the 
function  p(x,  y),  which  is  continuous  along  the  discontinuity  boundary  line. 
Two  possible  diffraction  wave  configurations  are  examined.  At  wedge  angles 
“1,  °2  <  a*>  the  limit  angle,  a  triple  configuration  (Fig.  1)  forms  near  the 


Fig.  1.  Diffraction  pattern  at  dj,  a  %  <  a 


intersections  L  and  F  of  the  shock  wave  with  the  edges  of  a  wedge.  The 
nonstationary  perturbation  from  the  shifted  wedge  center  O  propagates 
within  the  region  bounded  by  the  ABCD  arc  of  the  Mach  circle,  the  attached 
reflected  acoustic  wavefronts,  and  the  AD  portion  of  the  shock  wave.  At 
a  1 , 2  ®  Jr*  there  is  an  irregular  refraction  near  L  and  F  and  only  one  region 

of  nonstationary  flow,  which  includes  the  shock  wave  sections  outside  the 
wedge  (Fig.  2). 


fl r 


Fig.  2.  Diffraction  pattern  at  Otj,  a  2  >  a 
The  a*  value  in  the  system  of  x'y1  coordinates  is  determined  by 

ex*  —  arc  sin  (/.„  /  Y/ix2  -j-  Ay)  (h  —  Vi  —  A-) 


where  k0  =  D0/C0,  D0  =  shock  wave  velocity. 


(4) 


into  the  p.aneIr:yx;rr;/(Hr°3b,1rdldeitermining  P<X-^>  ls  transferred 

1  VI  \  ig.  and  the  luneis  represented  conformally 
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Fig.  3.  Region  of  non  stationary  flow 

in  the  Z  complex  plane. 


on  the  upper  half-plane  4  by  means  of  the  transform 

4  =  l  +  it]  =  +i») 


2«_1 


(5) 


The  analytical  function  F+,t,  of  t  half. piano  satisfies  the  condition 


where 


(6) 


1 

I 


o  =  0, 


«•  0 

S'. 


ft  -  tfl3 


#2  <KI<1) 

<m>i) 


f  !!!ti  Ih)^\ a  ~ ~  r*8  (6  ~ «*))  <V_  < «*) 


5i  —  In  -=  A,-1  (cosec  0,  —  /,• 
5s  —  \l~- 

and  8  (0)  is  a  Dirac  function. 


ctg°»).  It  —  lc  —  Ay1  (coscc  02  —  k  clg  02) 
h  ~  *  clS  «i.  5i  «  =  x-  ctg  a2 


(7) 


problem  with  continu^ paraLeters^  ^he^  bl^  !nh°m°eeneous  Hilbert 

the  function  having  a  zero  point  not  liwer  thZ^eT  8  'T**  **  determini"8 
solution  must  satisfy  the  condition  second  order  at  infinity.  The 


0  <*•(!- /><**,■  (i  +/) 


(8) 


Thtt  ^medium  ^  stationary 

function  medium.  The  boundary  <*>+(£)  of  the  analytical 


^+(a  = 


^2-/?2  +  iC  vrrp 


f 


s 


(9) 
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satisfies  the  condition  (8),  since  it  does  not  exhibit  singularities  in  the 
upper  half- plane.  The  function  F+({)/  $+  (£)  must  be  regular  at  infinity. 
It  follows  that 


(£)  =  O'  (?)  [V*  (£)  +  C0)  (10) 

where 


'I’MO  = 


1  1 

(  pt 

Pi 

~-T 

'1*T  I:,i)  (;* 

J 

//W  -1-/1,  Hi 

«!-»*+  A\lh' 

n 

\ •  *l-t 

s»-t  - 

)  ,  <*1.2  <(** 

Ch,2^>'*lt: 


(11) 


The  constant  CQ  in  (10)  is  determined  from  the  condition  of  smooth 
conjugation  of  the  wavefront. 

Pressure  is  determined  from 

p(6,ri)«Im$(^+i^)dt  +  /u  (12) 

-x 

where  p^  =  p^  for  a  j  ^  <  a  #  ar»d  Pa  =  ®  f°r  0  1  z  >  °  *• 

Pressure,  as  well  as  P ,  u,  and  w,  exhibit  a  logarithmic 
singularity  in  L  and  F.  A  formula  is  also  derived  for  p  distribution  p(y) 
along  the  curved  section  AD  of  the  shock  wave.  The  functions  u  and  w  are 
determined  from  (1)  and  expressed  in  terms  of  p  in  a  closed  system  of 
equations. 


It  is  noted  that  the  results  obtained  may  contribute  to  the 
development  of  a  theory  of  irregular  refraction. 


Salamandra,  G.  D.  and  N.  M.  Ventsel'. 

Measurement  of  gas  velocity  with  the  use  of 
the  contact  discontinuity  surface.  ZhTF, 
no.  11,  1971,  2463-2465, 

A  method  is  described  for  gas  velocity  measurements  behind 
a  detonation  wave  front.  The  contact  dis  continuity  surface  formed  in  the 
gas  is  used  as  the  point  of  reference  for  gas  motion;  the  contact  discontinuity 
surface  is  assumed  to  be  formed  by  the  interaction  between  shock  and 
detonation  waves.  The  surface  is  visualized  by  the  Toepler  technique.  For 
several  experiments,  the  numerical  value  of  gas  velocity  measured  16%  lower 
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^  ivsehrk  ,ube  \rrimeniB sh™  iha*  *«-  ««*«. 

v=iocily>  andy ; ziz^z::  \:ztz:zzz  i  shock  wave 

velocity  measurement  references,  e.  g.  the  heat  Inhum^e^' t^hn'iq™. 


Zaydel',  R.  M.  Shock  wave  passa^ 
tnrough  a  curved  boundary  interface 
between  two  media.  M7.Mr.  j  ]r}7-> 
111-121.  * 


a  light  medium  (X  > 

ii  during  the ' ■ 2  ^o/ou  (»,.  i,.  D\zx  izzv:zrPz? 


Fig.  1.  Shock  wave  interaction  with 
the  free- stream  interface  X  =  G  between 
two  media. 


the^oundary nd  “*  A  -  A’  «  at  a  time  X  along 


n  ~  <*.\p  (//,•  Y) 


(1) 


where  ka  Elisa  small  perturbation.  The  free-stream  boundary  is  set  in 
,X<  1  v'Jnd  ;  Va'reOC,hy  °  *°  *"  "eavy  Lldium 

shock  wale*,  A  set  of  H  '  Vel°citi”  °f  th«  “fleeted  and  transmitted 
ocr  waves.  A  set  of  gas  dynamic  equations  with  boundary  and  initial 
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conditions  is  formulated  for  small  perturbations.  Introducing  polar 
hyperbolic  variables 

y.  =r,chO„  x,  =  r.sh0.f  r,  =  y~y,2  —  x,\  th 6.  —  x,/y,  (2) 

and  using  a  Laplace  transform,  the  author  obtains  a  solution  to  the  equations 
in  the  general  form 

w"{q„  6.)  =  F.(q,  +  0.)  +  0.(g<  _  6.)  (s  =  1,  2)  (3) 

where  F  and$are  arbitrary  functions,  ps  =  sinhgs,  s  =  1  in  the  0<X  <V^ 
region,  and  s  =  2  in  the  -V2t  <X  <0  region.  Integration  of  the  original  gas 
dynamic  equation  in  polar  hyperbolic  coordinates,  using  (3),  gives  three 
finite  difference  equations  which  have  a  unique  solution  in  the  form  of  expansion 
in  an  exponential  series.  In  the  Ap — ►O  approximation,  where  Ap  is  the 
shock  wave  intensity,  the  solution  is 

^.(9.)  =  -}hv,°e-\  Ft(qi)  ==  '/2yt0e-«* 

(4) 

u(t)  =  katUty(i)  (x  a  icc,t) 

2 ni  J^r^)*pT^P==:  -f- 

9 

where 

V(p)  =p-Y,(p)  -f  ¥,(/>) 

V'(p)  =  (1  ~  p)  [p/M-  1/R]-',  W,(p,  =  (M  -  i/a)  [M  +  r //?]-. 

_  _  (6) 

V>(P)  —  (l-a)[l  +  pfl/r]-',  r=ypl-f-l,  fl=yay  +  1 

is  the  transform  of  'I'i(p),  C]  is  the  sound  velocity  in  the  s  =  1  region, P\tpj> 
and  a=  ci/c2.  Expression  {£>)  is  simplified  in  the  cases  of  a  =  1  (ci  =  c2). 

Pi  -Pi  (/*=  1)*  or  the  acoustic  impedance  pj,  c^  =  p., c2  (  Jt  =  a/U=  1). 
Integration  of  (5)  gives  the  asymptotic  valus 

oo)rs=  const  =  — — —  *=i  — — £i  (7) 

1  +  H  P.  +  P,  () 

i.  e.  'J'(t)  depends  on  the  P\l p-^  ratio  only.  A  generalization  of  the  conclusion 
led  to  an  equation  of  motion  for  the  free- stream  boundary 


<Pa 


d? 


P»  —  Pi 
Pi  +  Pi 


kga 


(8) 


where  the  acceleration  g  —  NU,  and  N  is  the  number  of  incident  shock  waves, 
each  imparting  AU  acceleration  to  the  boundary  interface.  This  equation 
is  identical  to  the  equation  of  Rayleigh- Taylor  for  gravitational  instability 
of  noncompressible  fluids.  It  is  concluded  that  the  effect  of  compressibility 
is  insignificant  when  the  wavelength  of  boundary  perturbations  is  sufficiently 
small. 


Shifrin,  E.  G.  Study  of  a  "hanging11  shock 
wave  near  the  point  of  origin.  MZhiG,  no.  6, 

1971,  30-37. 

Mapping  of  a  "hanging"  compression  shock  into  a  hodograph 
plane  is  described  for  a  two-dimensional,  nonuniform  supersonic  flow  in  a 
perfect  gas.  A  general  analytical  solution  of  the  shock  at  its  point  of  origin 
is  obtained  in  the  hodograph  plane  by  the  method  of  asymptotic  expansions  in 
first  and  second  order  approximations.  Shock  formation  conditions  are 
formulated  at  a  supersonic  point  of  flow.  It  is  shown  that  a  quadratic  parabolic 

form  of  the  boundary  line  in  the  hodograph  plane  is  a  necessary  condition 
of  shock  formation  at  a  supersonic  point  in  a  physical  plane  in  which  convolution 
is  present. 


Khristoforov,  B.  D.  Small-scale  electromagnet 
for  explosion  investigations.  PTE,  no.  1.  1972~ 

_  _*  i  '  9  '  9 


An  electromagnet  has  been  built  for  laboratory  measurements 
of  the  mass  and  wave  velocities  of  shock  and  detonation  waves.  The  field  of 
intensity  across  a  200  x  200  x  200  mm3  gap  is  ~  400  Oe;  in  addition  there 
is  a  permanent  magnetic  field  which  extends  over  a  50  x  50  x  40  mm3  region 

onnr  ^2^  Center*  The  70  mm  diameter  magnetic  circuit  is  made  of  900  x 
900  mm  Arraco  iron.  The  corners  of  the  device  are  strengthened  by  welding 
and  special  straps  to  increase  resistance  to  explosive  loads.  Dimensions  of 
the  magnet  gap  and  the  permanent  magnetic  field  are  controlled  by  specially 
designed  pole  shoes  made  of  a  200  mm  dia  type  St.  3  steel.  The  field 
intensity  across  the  gap  can  be  increased  in  inverse  ratio  to  the  end-area 
of  a  tapered  attachment  to  the  pole  shoe.  Two  paralleled  magnetic  coils 
each  composed  of  2000  turns  of  1.  8  mm  copper  wire,  are  used  for  superposed 
magnetization.  The  24  volt  current  supply  is  supplied  by  a  VSA6  rectifier. 
The  magnetic  circuit  is  installed  in  a  channel  welded  from  channel  bars*  the 
channel  itself  is  welded  to  a  table  made  of  angle  steel  to  increase  the  rigidity 
of  the  assembly.  Over  the  past  three  years,  the  electromagnet  has  withstood 
several  hundred  detonations  of  explosive  charges  up  to  0.  5  kg  without  sub¬ 
stantial  damage.  A  photograph  of  the  electromagnet  assembly  is  included 


Nevskiy,  L.  B.  Application  of  interferometer 
mirror  shift  for  gas  dynamic  investigations. 
OMP,  no.  2,  1972,  9-1  1.  - 


eas  Bow  intPrf^qo  1  1  ifntlySlS  lS  Preseuted  of  the  shift  of  supersonic 

grums’  Were  obtainefi  with  a  dual  beam  shift  inter- 

erometer  with  a  spherical  mirror  in  a  reflecced  divergent  optical  beam 
In  contrast  to  a  polarization  shift  interferometer,  this  interferometer  allows 
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a  smooth  and  uncomplicated  shift  of  the  wavefront 

□  CriVpH  fn-r  flrtiri  _  a  .1.  .  i  * 


Integral  equations  were 
incident  flow;  and  p  Ji  lu"  9  j°  ’  wt?ere  p  oo  is  the  density  of 

Experimental  data  obtained  from  the  shift 

WUH  the  Mach-V^er  Inte^™^  ^Xia^^nT  ^ 


Gorelov,  V.  A.,  and  L.  A.  Kil'dyushova.  An 
experimental  study  of  parameters  of  ionized 

17,  r  |rQnt  ?-'",tronR  shock  wave-  MZhiG,  no.  6, 


ahead  of  a  Stro!,g1h°c,ch  wav  "("f  “Vo  Ve'°City  “  i0"l2ed  air 

electric  discharge  shock  tube  «  t,  M,U1  pre."^  ‘n  “ 

w'tLX.T"'  X  f°r  *he  "«  ™*™«>  to  verify0 .iariier  dat'a  obtained"' 
by  the  authors  with  a  standard  probe  (MZhiG  no  2  1971  1471  a  , 

wire  shortcircuited  at  both  ends  was  used  Is  a  resonance  svstem  Tb  "lfCt°r 

smf  p  aCf'fl  al?n8  the  shock  tube  diameter  and  connected  to  the  feed  system” 
and  a  measuring ,  lute,  the  feed  system  was  connected  to  the  sh|!ge„era,„? 

V  ’  „  7  GH  )#  Resonance  Of  the  shorted  wire  appeared  at  the  instant 
tovZe^Th  WaS  1/2  N  X  ’  where  N  is  the  number  o  half-wave  on 

scope  "kce  The  C  r'C°rded  as  a  characteristic  spike  on  an  S" 

pe  trace.  The  corresponding  n  was  calculated  to  be  2.  8  x  1011  cm"3 

as  compared  to  the  2  2  *  lo"  cm. -5  vatue  measured  with  a  double  p^obe 
~  30V  „  I?  km/scc)-  The  tie  values  obtained  by  the  two  methods  differ  by 

wavefront  b",  "  “  =  .'<?>  P,ot’  »>>"«  *  is  the  distance  to  the  shock 

but  did°not  .ff«*  th™  rierveUJ'Cled  i"'°  th°  ",b<J  "ear  'he  P™* 

nri1„a„  ..  The  electron  mass  velocity  Ug  //  in  the  direction  of  shock  wave 
fntoPtahlatl°n  Wa%meaSUred  to  determine  the  effect  of  free -electron  diffusion 
was  <,  fe?l0n  Precursor  ionization.  An  electromagnetic  induction  method 
was  used,  based  on  measurement  of  the  potential  difference  at  the  boundary 
f  the  plasma  flow  through  a  transverse  magnetic  field.  The  experimental7 
He  versus  x  plots  show  that  the  measured  Ue  //  values  near  the  shoik  fron, 
hut %h  Sa,tlsfactorV  agreement  with  those  calculated  for  a  free  diffusion 
ut  they  decrease  rapidly  with  increases  in  distance  from  the  front  approach¬ 
ing  values  corresponding  to  ambipolar  diffusion.  PP 


j 

I 


1 


10 


Ivanov,  A.  A.  ,  L.  L.  Kozorovitskiy,  V,  D. 
Rusanov,  R.  Z.  Sagdeyev,  and  D.  P.  Sobolenko. 
j.xper imental  observation  of  electron  shock  waves 
in  a  collisionless  plasma^  ZhETF  P  v  TZ 
1971,  593-596.  ’  ' 


Experiments  are  described  which  establish  the  existence  of 
a  stationary  thermal  discontinuity,  or  electron  shockwave,  in  a  collision- 
less  plasma.  Tests  were  done  in  hydrogen,  argon,  and  xenon,  using  a 
p  asma  generated  in  a  glass  tube  by  two  shf  generators  and  an  axial 
magnetic  leld  in  the  0.  5--5  koe  range.  Local  heading  of  the  plasma  to 
electron  temperatures  of  some  300  ev  was  induced  by  a  single-turn  high 
current  cm1  generating  a  large  magnetoacoustic  wave  whose  energy  was 
absorbed  by  the  plasma  in  the  coil  region.  Probe  data  of  nT  vs.  axial 
position  then  show  a  drop  in  nT  and  pressure  characteristic  of  a  shock 
wave.  Variation  in  wavefront  parameters  were  investigated  under  different 
test  conditions;  these  showed  that  the  length  and  velocity  of  the  wavefront 
were  independent  of  initial  electron  temperature.  Tolerably  good  agreement 

was  found  between  theory  and  experiment  for  the  argon  and  hydrogen 

shock  data.  The  tests  thus  confirm  the  possibility  of  a  stationary  electron 
shock  wave  in  a  collisionless  plasma;  analogous  results  could  occur  from 
effective  braking  of  a  relativistic  high-current  beam  in  plasma. 
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II.  EQUATIONS  OF  STATE  FOR  GAS 


Stesik,  L.  N.  Calculation  of  detonation 
parameters  of  explosive  mixtures  with 
metals  using  an  equation  of  state  for  perfect 
gases.  FGiV,  no.  1,  1971,  111-117. 

Detonation  parameters  and  composition  of  detonation  products 
are  computed,  using  formulas  derived  for  a  system  satisfying  an  equation  of 
state  for  perfect  gases.  Detonation  velocity  D  is  determined  mainly  by  the 
relationship  of  two  opposing  factors:  the  heat  Q  released  in  a  detonation  wave 
and  the  fraction  of  a  condensed  phase  in  the  products.  Tabulated  theoretical 
data  and  plots  of  D  and  Q  versus  metal  percentages  led  to  the  conclusion  that 
the  D  of  the  metal  mixtures  (Al,  Be,  B)  in  explosives  having  a  negative  oxygen 
balance  (TEN,  pyroxylin,  RDX,  TNT)  increases  by  only  2-5%  with  an  increase 
in  metal  content,  while  D  in  mixtures  of  Al  with  oxydants  (hexanitroethane, 
ammonium  nitrate  or  perchlorate)  may  increase  by  10-30%. 


Ill .  EXPLODING  WIRES 


Yel'yashevich,  M.  a.  ,  A.  A.  Labuda,  L.  Ya. 
Min'ko,  I.  G.  Nekrashevich,  G.  M.  Novik,  and 
G.  I.  Bakanovich.  Pulsed  accelerator  generation 
of  hieh- velocity  plasma  beams  based  on  electric 


explosion  of  wires  and  dielectric  erosion. 


DAN  BSSR ,  v.  16,  no.  2,  1972,  115-117. 

An  improved  pulsed  generator  of  moving  plasma  is  described. 
Plasma  is  produced  by  an  electric  explosion  of  wires  in  a  discharge  chamber 
of  specific  geometry  and  dimensions.  In  one  version,  metallic  plasma  is 
produced  bv  electric  explosion.  Breakdown  of  current  flow  is  eliminated  by 
using  an  exploding  conductor  in  the  form  of- a  thin  metallic  film  deposited  on 
the  inner  surface  of  a  dielectric  insert  in  the  discharge  chamber.  Loss  of 
stored  capacitor  energy  is  avoided  by  placing  the  discharge  gap  inside  the 
chamber  (Fig.  1).  In  a  second  version,  a  pulsed  surface  dischaige,  without 


Ka  kv 


Fig.  1.  Pulsed  generator  and  typical  oscilloscope 
trace  of  a  discharge: 

1-  ring  electrode 

2-  dielectric  insert 

3-  central  electrode 

4-  firing  electrode 

5-  diagram  of  the  auxiliary  and  interna.1 
gaps  of  the  plasma  generator  using 
electric  explosion  of  wires 

netallic  film,  is  used  to  promote  erosion  plasma.  At  V  =  4kv  and  C  =  200mf 
n  the  first  version,  1.3  kj  was  discharged  for  a  1.6  kj  initial  energy.  _  The 
:alculated  energy  characteristics  of  the  plasma  generation  under  the  cited 
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conditions  are  plotted.  The  erosion  plasma  beams  produced  at  V  -  5kv  and 
C  =  1950  mf  were  recorded  by  high-speed  photography. 


In  both  versions  described,  underexpanded  hypersonic  plasma 
beams  can  be  produced  at  atmospheric  pressure.  A  low -temperature 
metallic  plasma  beam  of  a  predetermined  composition  is  proceed  wito  a 
relatively  high  efficiency  in  the  first  version  of  the  generator.  The  opt™ 
exit  velocity  of  plasma  was  25  km/sec.  A  relatively  dense  (ne  -  5  x  10  cm  ) 
low-temperature  (6000°  K)  erosion  plasma  is  formed  in  the  second  generator. 
The  erosion  plasma  exhibits  characteristic  high -intensity  continuous  and 

line  spectra. 
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HYPERSONIC  FLOW  PAST  BODIES 


Filatov,  Ye.  I.  The  optimum  shape  of  a  lifting 
body  at  hypersonic  speeds.  MZhiG,  no.  1,  1972, 

82-86. 

Optimization  of  the  aerodynamic  characteristics  of  a  horizontally 
flying  hypersonic  body  is  attempted,  omitting  some  previous  simplifying 
assumptions.  The  body  is  given  a  rectilinear  coordinate  system  in  which 
directions  of  the  x  and  z  axes  are  opposite  to  the  flight  direction  and  vertically 
downward,  respectively.  It  is  assumed  that  xz  is  the  symmetry  plane;  that 
z  =  my  and  z  -  f(x, y)  are  the  upper  and  lower  boundary  planes  of  the  semi-body 
y>0,  and  that  the  rectilinear  leading  edge  is  formed  by  intersection  of  the 
z  =  my  and  x  =  ky  planes.  It  is  further  assumed  that  the  bottom  section  is 
in  the  x  =  ki^=  i\  plane,  where  i7  is  the  semi-span;  the  function  f(x,y)  ••>n 
be  approximated  by  a  polynomial  of  a  low  order  in  any  x  =  const,  y  =  const 
section,  and  that  the  thickness  of  the  leading  edge  is  negligibly  small  in 
relation  to  the  characteristic  dimension  of  the  optimum  body  in  the  z  direction, 
i.  e.  at  x  =  ky  and  f(x,  y)  =  mz.  The  optimum  body  shape  problem  is  then 
formulated  in  terms  of  a- variational  method.  Solution  of  the  problem  is 
reduced  to  finding,  within  the  class  £  =  £  (£,tj)  of  surfaces  satisfying  the 
condition  C  (0,  tj  )  =  0,  a  surface  which  would  satisfy  the  maximum  value  of  the 
functional 

*=  (y,  +  y.)  (x,+xa)-‘  (1) 

where  Xj,  are  the  drag  factors  relative  to  the  double  velocity  head  on 
the  semi-body  and  its  leading  edge,  respectively,  and  Yj,  Y2  are  the 
corresponding  relative  lift  forces.  The  volume  of  the  semi-body  must  be 
V  =  const,  and  the  m,  !,  and  k  parameters  are  assumed  to  vary  with  surface 
geometry.  The  new  coordinates 


l  =  (*  —  ky)  1 1, 
nt= 

p  IF’ 


y/l,  £  =  (2  —  my) 1 1 
^L-pk  +  m 


are  introduced  in  the  solution.  Using  the  Ritz  method,  the  unknown  optimum 
surface  is  described  by  the  equation 


C“  8{6/*V+a..(Tl-0.25)  + 
+  812(11  —  0-1)  + 

+  020 (S  —  k  I  2)  -f  au  (|rj  — 

-&/l0)  +  a«(sTi*~ 
-k!  30)  +  ...} 
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The  problem  thus  formulated  contains  two  implicit  parameters:  the  coefficient 
of  friction  Cf  and  the  relative  radius  R  of  leading  edge  curvature.  Numerical 
calculation  of  the  unknown  coefficients  an,  a12,  etc.  in  (3)  was  made  by 
the  method  of  descent.  In  the  zero-th  order  approximation,  equation  (3) 
describes  pyramids.  For  a  pyramid  with  a  sharp  leading  edge,  calculations 
show  that  the  larger  the  span  (21),  the  higher  is  the  aerodynamic  quality. 

At  moderate  x-L,  the  quality  K  of  such  pyramids  and  the  slope  P  of  the  lower 
surface  toward  the  x  axis  only  weakly  depend  on  (Fig.  1  and  2).  At 
moderate  i\ ,  the  friction  T  is  one  third  of  the  total  drag. 


Fig.  1.  Maximized  parameters 
of  pyramids  with  sharp  leading 
edges  at  =  10-3. 


Fig.  2.  Maximized  para¬ 
meters  of  pyramids  with  sharp 
leading  edges  at  C{  =  10-2. 


.  ,  .  0  Trhe,  effect  of  bluntness  R  of  the  leading  edge  is  primarily 

evident  on  ^  of  the  optimized  pyramids  (Fig.  3).  The  parameter  P  remains 

constat  over  the  entire  range  of  R  in  contrast  to  the  pyramids  with  sharp 
leading  edges.  The  leading  edge  bluntness  does  not  significantly  affect 

the  aerodynamic  quality  of  a  lifting  body.  The  values  of  parameters  m,  /  k 

calculated  from  (3)  in  the  first  and  second,  approximation  indicate  that  the 
aerodynamic  quality  of  optimum  bodies  differs  by  less  than  0.  1%  from  that 
of  the  corresponding  pyramids.  The  dad  are  in  qualitative  and  quantitative 

JtklZrp  Teading'edges!1  We8tern  ^  appr°ximatinS  —-shaped  bodies 


.  *  18  concluded  that  the  geometry  of  the  body  type  studied, 

aving  a  given  volume  and  a  maximum  aerodynamic  quality,  is  nearly 
pyramidal  and  in  good  agreement  with  the  earlier  advanced  hypothesis 
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Fig.  3.  The  parameters  of  the  optimized 
pyramids  versus  blunt  R.  Continuous 
line  -  Cf  =  10“3;  broken  lines  -  Cf  =  10“2. 


Temkin,  L.  A.  Approximate  solution  to 
the  problem  of  the  ground  zone  in  a 
rarefied  gas.  MZhiG,  no.  1,  1972,  139-143. 

Two-dimensional  equilibrium  flow  in  a  monatomic  rarefied 
gas  around  a  semi-finite  rectangular  bar  is  analyzed,  on  the  basis  of  the 
Krook  model  for  collision  processes  in  gases.  The  Krook  equation  is 
given  in  the  form 


where  f  =  f(r,u)  is  the  coordinate  and  velocity  distribution  function  of  gas 
molecules,  fQ(r,u)  is  the  local  Maxwellian  distribution  function,  and  V 
is  the  integral  kinetic  operator.  A  set  of  equations 


*7.(r,u)du,  U(r)  =  -ij-JJJu.F/.(r,u)du 

T{t)  - 1 JTJ[u  - U(r)],VMr’  u>du 


(2) 
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was  derived  from  (1)  for  gas  density  n,  macroscopic  velocity  U,  and 
temperature  T.  The  set  of  equations  (2),  where  k  is  the  Boltzman  constant 
and  m  is  the  molecular  mass,  was  solved  by  iteration  in  a  first-order 
approximation  both  with  and  without  allowance  for  conditions  at  infinity. 

Using  the  M-20  computer,  the  macroscopic  flow  parameters  n,  U,  and  T 
wire  calculated  and  plotted  against  coordinates  x  and  y,  where  the  x-axis 
is  the  flow  direction  at  infinity.  The  plots  calculated  for  the  Knud  sen 
nvtnber  K  =  1  and  Mach  number  Moo  =  2,  and  for  K  oo  =  1  and  Moo  =  0.  75, 
displayed  a  sharp  rise  in  T  in  the  ground  zone  of  the  bar,  near  the  corners, 
followed  by  a  smooth  decrease  to  the  T  level  at  infinity.  T  on  the  axis  of 
symmetry  is  maximum  at  the  points  x  oo  ~  1  and  0.  6  at  M  ^  =2  and  0.  75, 
respectively.  At  both  Moo  values,  n  and  pressure  p=nT  increase  monotonically 
to  the  levels  at  infinity.  At  M  oo  =2,  the  stream  line  deviates  from  the 
body  near  the  ground  zone  because  of  decelerating  gas  flow.  Results  are 
in  qualitative  agreement  with  an  earlier  solution  to  the  ground  zone  problem 
obtained  within  the  framework  of  the  Navier-Stokes  equations. 


Shapiro,  Ye.  G.  Shock  layer  radiation  during 
hypersonic  air  flow  around  a  spherical  segment. 

MZhiG,  no.  1,  1972,  101-106. 

A  theoretical  analysis  is  presented  of  the  effect  of  shock 
layer  radiation  on  gas  dynamic  parameters  of  hypersonic  air  flow  (V »  = 
10-16  km/sec)  around  a  spherical  segment  of  a  base  v/ith  radius  R.  The 
shock  layer  flow  (Fig.  1)  is  described  by  a  set  of  integro-differential 


Fig.  1.  Shock  wave  profiles  and  flow 
line  system  for  a  segment  of  R  =  1  m  at 
V  oo  =  15  km/sec:  the  flow  line  numerals 
I  thru  V  correspond  to  0  =  0. 1 , 0.  3, 0.  5, 
0.7  and  0.9*  respectively. 
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equations  in  the  spherical  coordinate  system  r,  0  ,  whose  origin  is  shifted 
at  a  distance  /  =  0.  65  R  from  the  segment  center.  The  term  Q  or  the 
divergence  of  the  integral  radiation  flux  in  the  equations  was  apparently 
known.  The  equations  were  then  solved  for  gas  dynamic  parameters 
using  an  iterative  method  and  an  independent  variable 

l  —  (r  —  rb)  (r.  —  rb)~l  (1) 

where  rfi  is  the  shock  wave  profile  and  rb  is  the  body  profile.  All  gas 
dynamic  parameters  and  rs  were  approximated  by  Lagrangian  polynomials 
with  Hugoniot  relations  and  impermeability  as  boundary  conditions  in  the 
shock  wave  and  on  the  body  surface,  respectively. 

Calculations  were  made  for  segments  of  R  =  0-4  m  over  the 
cited  range  of  Vw  .  In  most  cases,  temperature  Too  ahead  of  the  shock 
wave  was  assumed  to  be  250°,  and  in  all  calculations  P,*  was  assumed  to 
be  10"4  atm.  As  with  flow  around  a  sphere,  the  radiation  effect  on  p  and 
V  was  found  to  be  practically  nonexistent  and  extremely  weak  with  respect 
to  the  distribution  of  the  velocity  component  v.  Also  analogous  to  the  flow 
around  a  sphere,  the  flow  difference  with  and  without  radiation  was  greatest 
at  the  body  surface.  The  temperature  distribution  parameter  was  the  one 
most  affected  by  shock  layer  radiation  (Fig.  2).  The  higher  the  flight  speed, 


Fig.  2.  Relative  temperature  distribution 
a£  Voo  =  9.85  km/sec,  R  =  1  m  with  (solid  lines) 
and  without  (dashed  lines)  allowance  for  radiation 
flow  along  th  lines  i>  2  =  0.  747,  i>4  =  1.141. 
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the  more  noticeable  was  the  radiation  effect  on  T  distribution.  Radiation 
cooling  was  also  intensified  by  an  increase  in  R  (Fig.  3)  The  plots  in 
Fig.  2  versus  those  in  Fig.  3  illustrate  the  effect  of  the  How  velocity  around 
the  segment.  A  decrease  in  T  of  the  shock  layer  results  in  a  decrease  in 


Fig.  3.  Relative  temperature  distribution 
at  V oo  =12  km/sec  and  R  =  0.  05-4. 

layer  thickness.  It  was  shown  that  the  radiation  energy  flux  to  the  stagnation 
point  increases  with  an  increase  in  R,  primarily  due  to  spectral  line 
radiation,  in  agreement  with  spectroscopic  theory.  The  universal  dependence 
of  the  ratio  of  radiation  flux  i;o  a  body  and  at  the  stagnation  point  is  also 
valid  for  a  spherical  segment.  An  earlier  conclusion  that  the  nose  of  a 
streamlined  body  is  one  area  where  radiation  heating  exceeds  convective 
heating  applies  also  to  a  spherical  segment. 

The  effect  of  radiation  near  the  rounding  point  on  gas  dynamic 
parameters  distribution  in  a  shock  layer  was  also  analyzed  on  the  basis  of 
the  energy  conservation  equation.  This  equation  indicates  that  as  Q  — ►  0,  the 
total  enthalpy  H  along  different  flow  lines  is  even  affected  by  small  radiation 
fluxes  (Fig.  4).  Fig.  4  illustrates  the  phenomenon  of  radiative  freezing 
which  occurs  when  a  particle  is  carried  out  of  the  region  affected  by  radiation 
to  a  colder  region  before  the  particle  energy  is  significantly  affected.  This 
phenomenon  determines  changes  in  distribution  of  gas  dynamic  parameters 
near  the  symmetry  axis  and,  in  particular,  causes  a  sharp  diop  of  H.  In 
regions  far  from  the  symmetry  axis,  particles  propagating  along  the  flow 
lines  do  not  lose  energy  because  of  the  high  velocity  and  low  Q  of  radiation 
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-  - ^uuun  ui  total  enthalpy 

H  along  different  flow  lines.  (  is  the 

distance  traveled  by  a  particle  from  the 
shock  wave. 


affect  distribution  of TeVodyna^ic^a ,  aP  ^  T  th®refore  does  not  significantly 
lations  showed  that  allowance  for  th  meters  of  the  shock  layer.  Calcu- 
not  required  in  the  vicinity  of  the  rtun^ng  poTnt^  radiative  volume  *■ 


Avduyevskty  V.  S. ,  V.  K.  Gretsov,  and 

ed'.esfan  *  with  leading 

^-^a11  zone-s-  M^hiC.  no.  1.  1972.  74-81. 

instability  waJrv^gTte^nb^  8ta11  zone 


Fig.  1.  a-  eemifinite  plate  with  fixed 

step  of  height  h  and  width  b;  b  -  cone  with 
a  fixed  shield  of  diameter  d. 
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isappearance  of  stall  zones  was  considered  as  instability.  The  geometry 
effects  of  the  streamline  bodies  on  flow  characteristics  were  studied  i  e 
the  parameters  /,  b/h,  and  p  of  the  plate  and  !  /d  and  0  of  the  cone! 

Shadow  photography  and  high  speed  motion  pictures  were  used  to  record 
gas  flow  phases.  An  unsteady  regime  with  pulsations  of  compression  shock 
was  observed  in  the  flow  past  the  plate  when  the  boundary  layer  in  front  of  the 
separation  point  was  either  laminar  or  turbulent.  The  pulsations  disappear 

at  b^h  p  =  70  or  at  !  <  A  ,  the  shock  wave  separation  value  in  a  perfect 
gas.  The  flow  stabilizes  when  b/h  is  decreased. 


The  flow  around  the  shielded  cone,  with  a  surface  laminar 
boundary  layer,  is  steady  at  a  sufficient  i>,  but  becomes  unsteady  when  j>/d 
is  decreased  below  a  certain  value.  The  flow  restabilizes  when  !  is  decreased 
o  a  /a  ue  <  ^.  Four  distinct  phases  of  pulsations  were  detected  in  the 

unsteady  flow.  The  flow  is  steady  at  0  >  0  ,•  ©  .  den<*nrifl  w  f 

Itrir,  ->\  rn  x  y  ^  ^  w  *  depends  on  M  of  free-stream 

flow  (Fig.  2).  Flow  stability  in  region  I  depends  on  i /d.  In  region  II,  the 


Fig.  2.  Angle  0  versus  Mach  number 
in  free-stream  flow.  Curve  1  -  theoretical 
®  *>  2-0  limits  of  separation. 


flow  is  steady  without  separation  of  the  turbulent  boundary  layer,  and  in 
region  III  the  shock  wave  is  separated.  The  value  of  the  angle  7  between 
the  Plate  or  cone  surfaces  and  a  line  connecting  the  leading  edge  or  cone 
tip  with  the  upper  edgc  of  the  step  or  shield  is  established  as  a  flow  stability 
criterion.  The  flow  is  steady  if  7  <  7*,  the  critical  value.  The  effect  of  V 

three-dimensional  flow  transition  on  pulsations  was  also  evaluated  using 
experimental  data.  K 


22 


Kokin,  G.  A,  and  Ye.  V.  Lysenko.  Pitot 
tube  pressure  in  supersonic  rarefied  gas 
flow.  MZhiG,  no.  1,  1972,  195-199. 

Readings  of  four  Pitot  tubes  in  a  low  density  wind  tunnel  at 
M  oo  =  2.  5,  2.6,  3.  2,  and  5.  4  are  compared  to  Pirani  and  McLeod  pressure 
gauge  data.  The  experiments  were  designed  to  evaluate  a  correction  to 
Pitot  tube  readings  in  a  supersonic  rarefied  gau  flow.  The  d/D  parameter 
of  the  experimental  Pitot  tubes  (Fig.  1)  was  0.  1-0.  2  and  0.  5-0.  9  for  1,  2  and 
3  tube  types,  and  the  X/d  parameters  was  1  and  50-100  for  the  1  and  2,  ^ 
types.  The  mean  free  path  \  oo  of  molecules  in  a  free  stream  flow  was  0.  5, 
1.2-0.  9,  3.  3-2.  4,  and  0.  5-0.  4  mm,  respectively,  at  the  cited  M  oo  values/ 

The  Knudsen  number  (Kn  oo  =  \  co  /D)  dependence  on  the  relative  Pitot  tube 
readings  at  the  cited  M  „  is  also  plotted  in  Fig.  1,  where  pOJ(c  is  the  measured 
Pitot  tube  pressure  and  po2  is  the  Pitot  pressure  behind  the  compression  shock. 


Fig.  1.  Geometry  of  Pitot  tubes  and  experi¬ 
mental  plots  of  P0*/P02  versus  Kn  oo  . 

Dotted  lines  -  theoretical  plots  for  r/R  =  0.5 
or  ~  0.863.  Arrows  indicate  theoretical 
values  corrected  for  thermal  conductivity  of 
the  tube. 

To  evaluate  the  correction  to  the  Pitot  tube  readings,  the 
correction  to  the  Pitot  tube  readings,  the  asymptotic  pressure  ratio  poF/po2 
in  free  molecular  flow  around  the  tube  was  formulated  using  known  relations, 
including  that  of  Rayleigh.  The  derived  formula  served  to  calculate  and 
plot  PoF/po2  ratios  for  thermally  insulated  tubes  of  different  l/d  parameters 
and  for  M  oo  =  1-9.  The  relation 
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was  derived  for  a  thermally  insulated,  long  cylindrical  tube  in  a  nearly 
free  molecular  streamline  flow,  where  S  «  =  M  K  /Voc  ,  the  ratio  between 
the  stream  mass  velocity  and  the  most  probable  thermal  velocity  of  molecules; 
and  r  and  R  are  the  inner  and  outer  radii  of  the  tube.  This  relation  was 
derived  in  approximation  of  a  single  collision  between  the  incident  molecules 
and  these  emitted  by  the  tube.  The  p0*/p02  values  for  Moo  =2.6  and  3.2 
calculated  from  (1)  are  in  good  agreement  with  corresponding  experimental 
values  (Fig.  1).  An  additional  correction  factor  C  for  thermal  conductivity 
was  introduced  into  (1)  to  obtain  theoretical  p0y^/p02  values  at  Moo  =5.4 
in  agreement  with  the  experimental  data.  The  experimental  and  theoretical 
data  were  used  to  plot  a  universal  Kn  dependence  of  the  parameter  £=  (p0*-p02) 
(PoF“Po2)  *  independent  of  M  oo  •  The  plot  can  then  be  used  to  determine 
Moo  from  measured  Pitot  tube  pressure. 
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Stulov,  V.  P.  jStrong  blowoff  on  a  blunt  bodv  surface 
in  supersonic  flow.  MZhiGT^'Z,  19 727^9^— 

flow  of  a  p0rfooethga°9rarou!,d  "aaM™t  bod™  wit°h  limuflaneou"^ 

surface  cliche injcc  °nK  »'  The  .wST'”  the  contact 

afth.  C  LrbodTe^eTThrf  (bOUndary  «■  «>e  .ho'ck  wave, 

spherical  a,8y8‘em  0l 

for  each  layer  assuming  that  the  pressure  on  iL  P  n  f  ^  ?°lvfd  seParat«ly 
the  Newton  formula  pressure  on  the  C  outer  side  is  given  by 


P it  —  Pi^sin5  o 


(1) 


where  <r  is  the  angle  of  the  flow  axle  with 
of  the  boundary  problem  by  the  method  of 
for  injection  according  to  the  formula 


C*  A  complete  numerical  solution 
successive  approximation  is  given 


U,  =  Uir0  cos’1  0 


(2) 


ov«hth!»e  fro^'aVrflct^n^^ra^c  ?0  TC“°d"  Ve,0blty 

cal  body  with  rc  radius  (Fig.  1).  The  discrepancy  bSw'een  the'approx 


Fig.  1.  Flow  diagram  around  a  blunt  body  B  with 
injection:  » 

1-  shock  layer,  2-  injected  gas  layer 


and  .ixact  numerical  solution  for  layer  2  is  shown  to  be  insignificant  Two 

™rirantTsabeasa"da  y“Cal  8°luti°n.  Variants  f°r  laV'  2  ««  Presented?  One 
l  and  2  la  b  .ed  ?n.t"  aBsumPll°"  of  constant  but  different  densities  of  the 
1?y'  *  In,thls  case*  two  equations  are  solved  for  rfi  and  r  .  The 
solution  tor  layer  1  approximates  well  the  numerical  data  calculated  ea- lie r 

neyar  or?a  7  n°W  ar°Und  d‘ffa""'  b°diaB-  Th'  -lotion  for Ifyerl 

solution  Th  T.l,  8y“me‘ry  “is  as,ees  wi,h  the  approximate  numerical 
(  _  ’  .  other  solution  variant  for  layer  2  was  obtained  within  the 

TV~  °*  b°™dary  layer  approximation  and  local  flow  self- similarity 
The  solution  in  this  case  agrees  with  the  first  variant  near  the  How  axis 
and  at  small  values  of  the  ratio  K  of  specific  kinetic  energy  of  iniection 
nd  incoming  flows.  But  the  discrepancy  with  numerical  solution  becomes 
significant  at  a  distance  from  the  axis,  e.g.  ~20%. 
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Mikhaylov,  V,  N.  ,  and  V.  S.  Tarmlov. 

Supersonic  flow  over  an  edge  formed  by 
intersecting  plates.  MZhiG,  no.  2,  1972, 

162  - 1 66. 

A  mathematical  method  is  introduced  for  the  numerical  snlnHnn 
of  a  supersonic  flow  problem  on  the  edge  formed  by  two  perpendicular  plates 
of  zero  thickness  In  a  Cartesian  coordinate  system,  the  plates  a^e  made 
to  coincide  with  the  y  =  0  and  z  =  0  planes,  and  the  velocity  vector  VcT^ 
defined  by  the  angle  of  attack  a  and  the  angle  0  between  the  v-axis  and 

Ind  ?heJvCfl°  >°f  V°°  °n  the  x  =  0  Plane  (Fig-  !)•  Pressure  p,  density  « 
and  the  velocity  components  u,  v,  w  must  satisfy  a  set  of  differentia/em.aHnnc 

with  plate  boundary  conditions.  Solution  of  the  equations  is  obtained  in  the 
region  bounded  by  the  y  =  1,  z  =  1,  and  x  >  0  planes  using  the  method  of 
adjustment  Formulas  are  Btve„  for  the  How  parameters  uv^wp  and 
p  and  the  stability  condition  of  the  calculation  scheme. 


Fig.  1.  Flow  diagram  over  the  edge  in 
perpendicular  plates 

,  _  ,MPlot.s  °f„the  calculated  parameters  along  the  half-lines 

Y;  .  •“  St/  co"^.al  variable)  show  that  the  parameters  become  constant 

a  perfect  'gas1*  presents  sh  A  ^ ualpative  theoretical  flow  pattern  plotted  for 
„,P  .  p  esents  shock  wave  traces  as  flow  region  boundaries  with 

sharply  different  parameters.  It  is  concluded  that  the  perfect  gas  model  is 
inadequate  to  describe  the  flow  type  studied,  because  of  a  discrepancy 
e  ween  the  pattern  and  an  earlier  experimental  pattern  of  interference  flow 
over  an  edge  in  intersecting  wedges.  The  method  of  calculation  was  also 
applied  to  the  expansion  flow  between  two  perpendicular  plates.  Calculations 

a  ,.ao  =  9  =  45°*  a  =  10°  show  that  flow  ahead  of  the  interference  region 

is  directed  from  the  edge  with  a  subsequent  significant  decrease  in  p  in  the 
interference  region. 


1.  Typical  plot  of  p  distribution  along 
the  y  axis  for  the  parameters  6, 

9  =  45°,  a  =  -10°  of  the  incoming  flow 
The  curves  1,  2,  3,  4,  5  correspond 
to  increasing  x  values 
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V.  ATMOSPHERIC  PHYSICS 


Bykov,  B.  P.  ,  and  S.  I.  Kozlov.  Ion  and 
electron  kinetics  In  a  disturbed  ionosphere 
at  altitudes  of  100^200  km.  Geomagnetizm 
i  aeronomiya,  no.  2,  1972,  340-342. 

Results  are  described  of  an  investigation  on  kinetics  of  ions 
and  electrons  at  altitudes  100  <  h  <  200  km  in  a  disturbed  atmosphere.  The 
disturbing  action  is  assumed  to  significantly  increase  the  degree  of  ionization 
a  of  air,  so  that  and  are  approximately  in  direct 

relation  to  the  N2,  O2  and  O  contents  at  these  altitudes.  Investigation  is 
limited  to  values  a  «  1,  so  it  is  considered  that  the  atmosphere  is  not  warmed 
up  due  to  a  sharp  rise  of  a^i.e.  electrons,  neutrons  and  ions  of  gas  have  the 
high  level  of  temperature  distribution  existing  in  normal  ionosphere. 


Form. 

no. 

Reaction 

Velocity  const  »cm^  per  sec 

1 

W'  +  Nj-  NO'  -|  N 

ai  =  2- 10'** 

2 

O'-I-O.  0.+  -1-0 

a-  =  4- !()-»» 

3 

NY  -|-  0  -►  Of  -f  Nj 

33  <  M0-“ 

4 

NY  +  0  —  NO1’  +  N 

si -- 2,5.10-'" 

5 

NY  +  0>  0/  +  Nj 

ctj  =  l.l  ()-w 

6 

OY  +  Nj  —  NO1'  +  NO 

a0  <  1  .  I0-'J 

7 

0  •+  +  «  -  0  -!-  0 

a*  (0»+)  =  2,2.  in-7 (7’e/300)-i 

8 

NY  -1- «  -  N  +  N 

a*  (NY)  =  3- 10-7  (/’f/3(l0)‘,/1 

9 

NO*  +  «  N  -1-  0 

a*  (NO ')  =  4 ,4  •  10-7  (7,</300)-»,  200"  K  < 

<  7’e  <  500’  K 

a'(NO')  =  0,9- 10-7  (7’c/1000)-*,«,  Te  >  500°  K 

Table  1. 


The  general  scheme  of  elementary  processes  used  in  this  work 
and  their  velocity  constants  are  given  in  Table  1.  On  the  basis  of  these 
reactions  kinetic  differential  equations  are  obtained  which  were  solved  by 
computer.  Calculations  are  done  for  a  wide  range  of  altitudes  and  values  of 
Neo»  day  and  night  conditions.  Results  of  calculations  for  daylight 
conditions  at  <13  =  10“  ^  and  <X£,  =  10“  ^  cm^/sec  are  plotted.  The  main 
characteristic  features  of  ion  and  electron  kinetics  in  a  disturbed  ionosphere 
are  seen  to  be  as  follows:  1)  The  role  of  ions  is  significant  only  at  the 
initial  moment  after  disturbance,  not  exceeding  in  most  cases  —1  sec; 
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2)  At  lower  altitudes,  as  in  a  normal  atmosphere,  recombination  takes 
place  mainly  according  to  a  square  law.  However,  with  an  increase  of 
Neo  after  the  disappearance  of  N2  over  a  certain  time,  velocity  variation 
may  follow  a  linear  law.  The  duration  of  this  time  significantly  increases 
at  higher  altitudes;  3)  The  time  behavior  of  [NO+J  differs  little  from 
predictions.  Even  at  upper  altitude  regions,  the  role  of  this  ion  is  small, 
except  when  Neo  <  10  cm-3.  4)  To  a  first  approximation,  at  all  altitudes 

the  main  molecular  ion  is  Oj,  ;  5)  The  principal  difference  in  behavior  of 
10ns  and  electrons  in  night  conditions  from  day  consists  in  the  fact  that 
the  dissociative  recombination  level  increases,  so  that  the  role  of  N  +  and 
NO  ion.  substantially  decreases;  and  6)  The  decay  time  of  electron  concen 
tration  from  Neo  to  normal  values  of  an  undisturbed  ionosphere  has  only  a 
slight  dependence  on  the  value  of  Neo. 

These  conclusions  also  hold  up  under  varying  values  of  the 
unknown  velocity  constants  a 3  and  a6>  All  calculations  were  done  assuming 
U1  -  to  be  independent  of  temperature. 


Alimov,  V.  A.,  L.  M.  Yerukhimov,  and  T.  S. 

Pyrkova.  Theory  of  the  F gprpafj  phenomenon 
in  the  ionosphere.  Geomagnetizm  i  aeronomiya, 
no.  5,  1971,  790-797. 

An  analysis  is  given  of  possible  factors  which  cause  stretching 
of  the  return  r-f  pulse  from  an  ionospheric  probe,  or  the  F8Dread  effect. 

The  authors  use  the  diffraction  theory  of  pulsed  signal  spread  in  a  nonuniform 
ionosphere,  together  with  data  on  nonuniformity  distribution  vs.  altitude 
in  the  F-layer,  and  diffusivity  data  on  returned  signals.  Both  ground-based 
and  satellite  transmitters  are  considered,  and  the  variation  in  spread  effect 
at  different  latitudes  is  considered.  Using  an  averaged  model  of  the  F-layer 
structure,  the  authors  obtain  expressions  for  electron  density  variation 
AN/N  in  the  equatorial,  middle  and  high  latitudes.  Following  this  the  character¬ 
istic  of  returned  pulsewidth  T  is  obtained  for  several  possible  AN/N  models, 
assuming  either  an  exponential  or  parabolic  variation  within  a  layer,  as 
tabulated  in  Table  1.  Graphical  comparisons  are  given  of  theoretical  and 
calculated  behavior  of  the  depth  of  the  diffusivity  zone  Az  =  c  t  as  a  function  of 
signal  frequency.  Fig.  1  gives  this  relation  for  the  extremes  of  equatorial 
and  polar  latitudes,  based  on  ionogram  data  from  the  Alouette  satellite.  These 
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Model  of  non- 
uniform  laver 

Exponential 

layer 

Parabolic 

layer 

"(A7V)*(=)  =  const 

T  ~  O'* 

T  ~  (0 

(&NIN)1  (*)  =  const 

r  const 

T~0>* 

(5/Vj*  (i)  =  D/A  —  t 

— 

T  m  const 

(UN /N)s  (*)  =  Aie~al 

T  ~  0) 

— 

(AJV//V)’(*)  =  Ci  + A«az  ■ 

T  ~  V  Oi  +  ih/(i>’ 

Table  1. 


Fig.  1.  Az(f)for  satellite  signals, 
a  -  equatorial,  b  -  polar  latitudes 

results  show  a  distinct  variance  from  analogous  data  using  a  ground-based 
probe  system,  as  seen  in  Fig.  2,  indicating  the  discrete  nature  of  inhomo¬ 
geneities  in  differing  portions  of  the  F-layer.  The  findings  appear  to 


Fig.  2.  A  z(f)  for  ground-based  signals, 
a  -  equatorial,  b  -  polar  latitudes 


<? 
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establish  that  the  Fspread  mechanism  is  generally  caused  by  scattering 
o  he  signal  pulse  on  local  electron  den.it,  variation,  in  area,  on  "e  order 
Of  l-5km  across.  The  variation  in  spread  behavior  is  then  governed  bv 

°f  'he  “  *•  reflected  SL 


Alimov,  V.  A.  Frequency  correlation  of  fluctuations 
in  radiowaves  reflected  from  the  ionosphere.  Geo^ 
magnetizm  i  aeronomiya,  no.  3,  1972,  548-551. 

.  ,  Jhe  author  considers  the  question  of  frequency  correlation  in 

the  received  fluctuations  of  dual  frequency  waves  reflected  from  the  iono¬ 
sphere,  as  a  means  of  interpreting  some  inhomogeneous  characteristics  of 
the  ionosphere.  As  noted  previously,  the  correlation  function  R  (f  f  *) 
may  vary  both  within  the  ionosphere  and  beyond  it;  also  vertical  W1  w2 
r?fra^lonJ  effects  must  be  allowed  for.  Alimov  analyzes  the  general  case 
of  inclined  incidence  of  the  transmitted  beam,  and  determines  expressions 
r  l-  *“2**  at  varYing  conditions  of  ionospheric  anisotropy.  Analysis 
of  this  case  generally  requires  calculation  of  the  trajectories  of  wave  normals 
when  a  substantial  effect  of  vertical  refraction  on  R  (f^f^)  is  to  be  expected. 


Vul  fson,  N.  I.  ,  and  L.  M.  Levin.  Explosive  breakup 

-developing  cumulus  clouds.  FAiO~o'._2,  1972, - 

156-lfa6. 


Upward  and  downward  spontaneous  convective  jets  initiated  bv 
explosKms  at  various  heights  in  developing  cumulus  clouds  were  investigated. 
Re  ationships  are  found  among  the  parameters  of  jets  formed  in  different 
sectors  of  an  unstable  layer  (a  conventional  developing  cumulus)  to  determine 

Z^neS  m  ah  exploslons  create  significantly  more  intensive  downward 
rather  than  upward  movements.  Mathematical  expressions  are  derived  and 
results  are  tabulated.  Comparison  of  the  :alculaied  jet  velocities  shows 

of^nnntl0810118  w  ^  uPfer'  sectors  developing  cumulus  produce  a  system 
of  spontaneous  jets  with  a  destructive  capacity  (caused  by  the  downward  jets) 
considerably  exceeding  the  intensity  of  cloud  development,  due  to  upward  jets 
The  more  favorable  the  condition  of  cloud  development,  the  greater  the 
intensity  of  cloud  destruction.  Experiments  were  conducted  under  natural 
conditions  in  the  I  ergansk  valley  during  May  1970,  using  explosive  shells 
in  dense  convective  clouds.  The  shooting  of  two  cumulus,  5500  and  6600  m 
hick  by  antiaircraft  mine  shells  lowered  their  thickness  by  3  to  4  times  and 

a  u  gradually  va^lshed-  The  explosion  effectiveness  can  be  further  in¬ 
creased  by  blasting  with  a  special  type  of  cumulative  shell,  which  generates 
downward  jets  during  explosion.  generates 
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VI.  SOIL  MECHANICS 


Iskakbayev,- A.  Crack  propagation  in 
the  curve  of  a  linear  viscoelastic  fold. 

VAN  KazSSR,  no.  12,  1971,  54-57. 

Brittle  fracture  of  a  viscoelastic  layer  (rectangular  band) 
compressed  at  the  ends  by  a  force  p  within  a  fold  of  a  rock  mass  is  analyzed. 
The  fracture  is  examined  as  a  jointing  process  with  a  creep  deformation 
background.  The  material  of  the  rock  mass  is  described  by  the  standard 
linear  body  model.  The  time  required  for'  fracture  initiation  is  determined 
from  a  formula  of  tensile  stress  in  the  curved  section  of  a  fold  subjected  to 
a  bending  force  (p-q),  where  q  is  the  weight  of  the  rock  mass.  Solution  of  a 
transcendental  equation  gives  the  time  t  >  tj  necessary  for  fracture  propagation 
through  a  28  thick  layer  cross  section  from  a  given  point  of  the  surface  fiber. 
Determination  of  the  time  of  layer  fracture  is  similarly  arrived  at  when  the 
layers  are  described  by  a  Maxwellian  model. 


Rodionov,  V.  N.  ,  V.  V  4.dushkin,  V.  N. 

Kostyuchenko,  V.  N.  N...olayevskiy,  A.  N. 

Romashov,  and  V.  M.  Tsvetkov.  The  mech¬ 
anical  effect  of  an  underground  explosion. 

Moskva,  Izd-vo  Nedra,  1971,  224  p. 

A  comprehensive  review  is  given  of  all  fundamental  character¬ 
istics  relating  to  mechanical  effects  of  underground  explosions.  A  simplified 
method  is  set  forth  for  calculating  explosion-produced  mechanical  motion , 
which  enables  prediction  of  explosion  cavity  dimensions,  destruction  range, 
fissurability  of  rock,  pa  rameters  of  cratering  explosions,  and  intensity  of 
seismic  waves.  Theoretical  results  are  compared  with  experimental  labora¬ 
tory  and  field  data  from  Soviet  and  non-Soviet  sources.  Engineering  applica¬ 
tions  of  strong  underground  explosions  are  discussed;  examples  given  include 
the  creep  dam  against  mud-debris  flows  formed  in  Medeo,  near  Alma-Ata, 
in  19b6-b7,  and  the  Baypazinsk  explosion  of  1968  which  generated  a  rock 
slide  to  dam  the  Vakhsh  River.  Theoretical  solutions  to  analogous  problems 
are  included,  using  basic  physical  properties  of  soil  and  rock,  and  possible 
future  developments  in  the  techniques  are  suggested. 
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VII.  LASER  SIMULATION  AND  RELATED  EFFECTS 

\ 

Ashmarin,  I.  I.,  Yu.  A.  Bykovskiy,  V.  A. 

Gridin,  V.  F.  Yelesin,  A.  I.  Larkin,  and 
I.  P.  Sipaylo.  Shock  waves  generated  by  the 
action  of  laser  radiation  on  transparent  bodies. 

IN:  Sbornik.  Kvantovaya  elekti  onika.  Moskva, 

Izd-vo  Sovetskoye  radio,  no.  6,  1971,  126-128. 


A  set  of  experiments  is  described  on  interaction  of  powerful 
laser  radiation  with  type  X-8  glass,  with  the  object  of  determining  the 
criteria  for  shock  wave  generation  at  the  focal  point.  A  Q- switched  ruby 
was  used  generating  50  Mw  pulses,  focused  to  a  0.2  mm  spot  both  internally 
and  on  the  glass  surface.  A  holographic  plus  high-speed  framing  method 
was  used  to  determine  wave  propagation  velocity.  Results  confirmed  that 
the  transition  from  a  longitudinal  sonic  wave  to  a  shock  wave  generally 
occurred  near  the  damage  threshold.  Wave  velocity  was  observed  to  be  a 
function  of  pulse  intensity,  as  well  as  of  delay  time  in  pulse  application, 
as  seen  in  Fig.  1.  In  theory  the  sonic  wave  converts  to  a  shock  wave  when 


I,rel.  units  \ 


0,3 

0,6 

0, * 
0,1 
0 


A£,n  S 


i-o-l 

KH 


U3-1 


160 


*120 


SO 


K>l 

_ L_ 


7  0 
V,  10  s  cm/  sec 


Fig.  1.  Wave  velocity  vs.  delay  (squares) 
and  laser  pulse  intensity  (circles). 


pressure  p  in  the  focal  region  attains  a  value  of  pQ Cq  ,  where  pQ  is  glass 
density  and  CQ  is  the  sonic  velocity  in  the  glass.  A  sample  calculation  based 
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on  actual  parameters  was  made  on  this  assumption,  showing  that  shock 
wave  velocity  should  be  1.3  times  sonic;  this  agreed  closely  with  measured 
results.  With  the  beam  focused  on  the  surface  a  dual  spherical  wave  was 
typically  generated,  as  shown  in  Fig.  2.  Their  characteristics  were  essentially 
the  same  as  for  the  internally  focused  pulse. 


Fig.  2.  Dual  wave  generation  from 
surface  irradiation. 


Cherkun,  Yu.  P.  ,  and  I.  N.  Konopel'ko. 

Current  status  of  experimental  and  theoretical 

studies  in  the  field  of  combined  heat  exchange. 

I-FZh,  v.  22,  no.  4,  1972,  757-758. 

A  brief  summary  is  given  of  several  papers  on  the  subject 
of  radiation  -  induced  heat  exchange,  presented  at  a  session  of  a  general 
conference  on  mass  and  heat  transfer  held  on  Sept.  15-17,  1971  at  the 
Thermophysics  Institute,  Siberian  Branch  AN  SSSR.  The  general  theme  of 
the  session  was  phenomena  related  to  combination  heat  transfer  mechanisms, 
namely  radiative  -  conductive  and  radiative  -  convective. 

S.  S.  Kutateladze  and  D.  I.  Avaliani  reported  on  experimental 
findings  of  laser  beam  attenuation  in  liquid  media  with  pulsating  turbulence; 
this  work  was  judged  important  from  a  heat  transfer  viewpoint  as  well  as  in 
regard  to  beam  propagation  in  a  turbulent  medium.  A  related  paper  by 
V.  M.  Kostylev  and  V.  Ya.  Belostotskaya  dealt  with  the  relation  between 
conductive  and  radiative  heat  transfer  in  optically  thin  layers  of  finely-dispersed 
media,  they  shew  that  in  this  case  it  is  incorrect  to  use  the  hypothesis  of 
additive  fluxes. 

A  report  by  N.  A.  Rubtsova  and  A.  E.  Berte  covered  the 
thermal  state  of  pure  metal  surfaces  exposed  to  a  constant  thermal  flux. 

Results  showed  anomalies  in  surface  temperature  distribution  (the  Jacques 


Reproduced  from 
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effect)  on  metals  including  Al,  Fe,  Sn  and  others; 
for  this  phenomenon  are  advanced. 


several  possible  models 


....•i  V-  M*  Yeroshenko  discussed  variants  of  ablation  techniques 
to  this  Zf  °/  dlSPCrSCd  mCdia  iS  UScd  as  a  protection  against  radiative 
20  no  Sc  l  9VOCesSeS  at  level*  of  100  atm  and  temperatures  in  the 

20,000  C  range.  Several  papers  by  other  authors  dealt  with  heat  exchange 
studies  in  various  types  of  combustion  reactions.  xcnange 

A  M  p  i  Sf?CCial  mention  is  given  to  a  paper  by  A.  I.  Leont'yev  and 
A.  M  Pavlyuchenko  which  is  cited  as  giving  the  most  complete  picture 
to  date  on  temperature  distribution  in  a  turbulent  boundary  layer  and  laminar 
sub.ayer  »  a  thermal  radiation  field.  As  with  Kostylev  eLlJ  tlZ  ,Z 

clS„dit^“Stra‘e  *he  faUaCy  °f  ‘he  additive-fl“x  hypothesis  under  their  test 


in  closing,  b.  Kutateladze  stressed  the  need  for  further  exn 
mental  data  on  radiative  heat  exchange  resulting  from  interaction  of  e-m~' 
waves  with  materials.  e  m 


Epshteyn,  E.  M.  Thermal  instability  of  a 
semiconductor  in  a  laser  beam.  IVUZ  Radiofiz 
no.  1,  1972,  33-37. 


i  •  ^  A  theoretical  study  is  presented  on  a  possible  mechanism  of 
laser- induced  destruction  of  a  semiconductor.  The  case  is  limited  to  photon 
nergies  less  than  the  width  of  the  forbidden  zone,  so  that  the  main  exchange 
echanism  is  absorption  by  free  carriers.  In  the  general  case  a  cumulative 
heating  and  increase  in  free  carriers  will  occur,  resulting  in  a  stationary 

f.  "\P®ra  ,re  ,field*  ^OWever’  above  some  laser  threshold  level  the  stationary 
leld  breaks  down  and  a  rapid  thermal  reaction  sets  in  which  is  analogous  to  > 
breakdown  in  a  d-c  field,  or  the  thermal  explosion  from  an  exothermic 
chemical  reaction.  To  examine  the  governing  phenomena  the  author  assumes 
oimplest  case  of  a  cylindrical  semiconductor  target,  axially  excited  by 
a  uniform- intensity  laser  beam  such  that  l«  p  «  R>  where  X  is  laser 
wavelength,  and  p  and  R  are  beam  and  cylinder  radii.  Expressions  are  derive 
for  the  limit  conaitions  of  the  stationary  field  solution,  and  the  corresponding 
threshold  power  of  the  laser  is  determined.  Finally,  it  is  shown  that  the 
elapsed  time  required  to  reach  threshold  under  practical  conditions  will  lie 
in  e  millisec  -  1  sec  range  -  i.e.  demanding  in  effect  a  c-w  laser  regime 
Therefore  the  described  mechanism  cannot  account  for  breakdown  observed 
for  example,  in  transparent  dielectrics  under  nanosecond  pulse  exposure  ' 
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Frolov,  V.  V.  Temperature  fields  in  multi¬ 
layer  translucent  coatings  under  conditions  of 
pulsed  radiation  heating.  I-FZh,  v.  22,  no.  4, 

1972,  755-756.  {Annotation  of  deposited  paper.  ) 

A  summary  is  given  of  an  analysis  of  heat  transfer  in  a  multi¬ 
layer  coating  system  under  beamed  irradiation.  The  problem  is  treated 
as  nonstationary  and  one-dimensional  with  distribution  of  absorption  taken 
into  account  Reflection  and  concentrated  absorption  at  layer  interfaces 
are  treated;  the  absorption  process  is  considered  to  be  linear. 

As  an  example  the  temperature  distribution  is  calculated  in 
a  two-layer  heat  shield  in  which  the  heat  load  is  assumed  to  be  applied  in 
triangular  pulses  of  a  given  energy,  at  durations  ranging  from  10  to  60 
seconds.  Owing  to  changing  boundary  conditions,  it  is  shown  that  the  system 
must  be  solved  as  a  class  of  piecewise-discontinuous  functions.  On  this 
basis  graphical  solutions  are  obtained  showing  the  time  characteristic  of 
heat  absorption  as  a  function  of  absorption  coefficient  and  pulse  duration. 
Instantaneous  temperature  distributions  are  also  presented  for  the  cases  of 
absolutely  opaque  as  well  as  semitransparent  coating  systems. 


The  analysis  shows  that  failure  to  account  properly  for  internal 
heat  absorption  can  lead  to  significant  errors  in  theoretical  temperature 
fields  and  level  of  total  absorption;  furthermore,  the  error  increases  with 
reduction  in  pulse  width. 


Bessarab,  Ya.  Ya.  ,  Yu.  B.  Tkach,  V.  P. 
Zeydlits,  N.  P.  Gadetskiy,  and  V.  V.  Dyatlova. 
Study  of  collective  processes  in  a  plasma  using 
light  of  stimulated  emission.  Z!  £TF,  v.  62, 
no.  2,  1972,  569-572. 


The  use  of  coherent  light  emission  to  study  the  collective 
processes  in  a  plasma-beam  discharge  is  described.  It  is  shown  that  an 
e  ective  method  of  studying  these  processes  is  the  direct  observation  of  the 
orm  of  oscillation  in  the  intensity  of  spontaneous-emission  spectral  lines, 
with  a  subsequent  application  of  Fourier  analysis.  Study  of  the  oscillation 
spectrum  was  performed  through  investigation  of  the  fluctuation  intensity  of 
coherent  emission  in  the  Ar  II  transitions  4p2D°5/2  — *  4S2p3/2  at  a  wavelengt 
f  4880  A,  and  4p*D%  „  -  4s2p3/2  a,  a  wave^th  of  51«*£  The 
frequency  excited  in  the  plasma,  which  depends  upon  modulation  of  the 
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emission  intensity,  is  determined  by  the  relationship  o>maxT^l  where  T  is 
the  lifetime  at  the  upper  level  of  the  observed  transition.  In  this  study 
(i)  max  ~  100-120  MHz.  In  the  case  of  coherent  emission,  the  spectrum 
of  the  oscillation  being  investigated  has  a  lower  limit  given  by 


0)njin 


2zE2c{i—R) 
(Ni — iV:)/io)ud  ’ 


where  eQ  is  the  dielectric  permeability  of  the  active  medium,  E  is  the  field 
intensity  of  the  light  wave  in  the  resonator,  c  is  the  speed  of  light,  R  is  the 
reflection  coefficient  of  the  resonator  mirrors,  d  is  the  distance  between 
mirrors,  1;  is  Planck's  constant,  is  the  frequency  of  the  observed  transition 

and  and  N2  the  corresponding  population  densities  of  the  upper  and  lower 
laser  levels,  respectively.  Fig.  1  compares  obtained  spectra  of  coherent 
and  spontaneous  emission. 


S,  rel.  units 


Fig.  1.  Realization  and  type  of' spectrum 
obtained,  a  -  by  modulation  of  coherent 
emission,  b  -  by  modulation  of  spontaneous 
emission. 


Sultanov,  M.  A.  Destruction  of  transparent 
dielectrics  under  the  action  of  free-running 
neodymium  and  ruby  lasers.  Mekhanika 
polimerov,  no.  2,  1972,  359-360. 

The  author  notes  that  a  variety  of  theories  exist  on  the  basic 
destruction  mechanism  of  laser  interaction  with  transparent  soiids;  among 
these  are  thermal  explosion,  acoustic  phonons  from  stimulated  Brillouin 
scattering,  and  photoiy  ci®  of  the  target  material.  Experiments  are  described 
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which  were  designed  to  clarify  the  predominant  mechanism  of  laser  damage 
in  several  dielectrics,  including  polymethyl-methacrylate,  polystyrene  and 
several  types  of  glass.  Both  ruby  and  neodymium  crystal  lasers  were  used, 
in  free-running  regimes  at  energies  of  3  to  15  j.  Characteristic  platelets  of 
cracks  at  angles  up  to  90c  from  the  beam  axis  were  observed  for  both  laser 
types;  however  in  the  Nd  target  specimens  the  damage  area  was  more 
extended  along  the  beam  axis,  whereas  the  ruby  damage  was  localized  near 
the  beam  focus.  (Repeated  reference  is  made  to  damage  photos,  which  are 
unfortunately  omitted  from  the  text.  ) 

The  test  results  lead  Sultanov  to  the  following  concept  of  the 
damage  p.  jcess.  On  focusing  a  sufficiently  powerful  beam  in  the  transparent 
dielectric,  a  local  region  is  heated  to  a  liquid  state  in  which  further  energy 
is  absorbed  by  inhomogeneities,  and  crack  propagation  begins.  A  plasma  is 
generated  which  expands  at  high  velocity;  thus  the  effect  is  that  of  a  powerful, 
instantaneous  point  explosion.  The  process  is  a  hydrodynamic  one,  including 
shock  wave  generation  which  causes  further  destruction  both  before  and 
behind  the  beam  focal  point.  Examples  of  remote  damage  regimes  in  PMMA 
specimens  are  cited  in  added  support  of  the  theoretical  model. 


Korotin,  A.  V. ,  and  L.  P.  Semenov.  Vaporization 
of  crystals  under  the  effect  of  external  excitation. 

IN;  Institut  eksperimental'noy  meteorologii.  Trudy, 
vyp.  30.  Fizika  aerodispersnykh  sistem.  Moskovskoye 
otdeleniye  gidrometeoizdata,  Moskva,  1972,  65-71. 

The  authors  present  a  straightforward  thermodynamic  analysis 
of  the  interaction  of  a  concentrated  heat  flux  with  a  crystal  surface,  for  the 
case  in  which  an  appreciable  melt  zone  appears  prior  to  evaporation.  The 
analysis  assumes  a  constant-intensity  beam  normal  to  a  semiinfinite  crystal 
face,  and  arrives  at  expressions  for  melt  zone  boundaries,  growth  rate 
and  limit  conditions,  and  time  to  melt,  in  terms  of  target  and  beam  para¬ 
meters.  It  is  shown  that  in  the  general  case  the  maximum  temperature  will 
occur  at  the  outer  melt  surface,  and  also  that  an  optimum  beam  intensity 
exists  for  which  the  melt  area  will  be  maximum,  decreasing  at  higher  or 
lower  intensities.  It  is  interesting  to  note  that  the  numerical  examples  given 
assume  an  ice  target;  results  on  tie  Ice-water  parameters  are  given  for  beam 
densities  ranging  from  25  to  200  w/cm  ,  In  ice,  for  example,  the  depth  of 
the  melt  zone  is  only  weakly  dependent  on  beam  intensity. 
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Basov,  N.  G.,  Yu.  S.  Ivanov,  O.  N.  Krokhin, 

A.  Mikhaylov,  G.  V.  Sklizkov,  and  S.  I. 
Fedotov.  Generation  of  neutrons  from  spherical 
irradiation  of  a  target  by  powerful  laser  radia- 
tion»  ZhETF  P,  v.  15,  no.  10,  1972,  589-591, 


*  u  ^  T^e  authorS  note  some  limitations  to  neutron  production  from 
aser  heating  of  a  target  for  fusion  purposes.  Specifically,  the  effect  of  an 
increasingly  powerful  focused  laser  becomes  offset  by  diffusion  of  the  high 
temperature  region  owing  to  thermoconductive  and  gasdynamic  enerev  loss 
An  alternative  approach  suggested  recently  by  Basov  et  al  is  to  heat^a 
spherxcal  target  simultaneously  with  multiple  beams;  in  the  present  case 
this  was  done  with  a  deuterated  polyethylene  target  exposed  to  nine  equal 
beams,  as  indicated  in  Fig.  1,  using  an  Nrf  glass  laser  in  the  giant  pulse 
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Fig.  1.  Multibeam  array  for  CTR  target. 

1-  preamplified  beam 
6-8-  second  amplifier 
21-29-  third  amplifier 
42-50-  focus  lenses 

Compensating  delays  for  differing  path  lengths 
not  shown. 


mode.  This  array  attained  a  mean  power  density  of  1016  w/cm2  on  the  target 
surface,  at  2--16  ns  duration.  The  focusing  objectives  were  placed  to  obtain 
a  focal  plane  200g.  from  the  target,  for  minimum  reflection  and  uniform 
heating. 


Some  results  are  shown  in  Table  I  for  various  target  sizes  and 
beam  energies;  the  measured  value  was  obtained  from  three  scintillation 
counters.  The  nr  values,  calculated  independently  for  thermoconductive 
and  gas  dynamic  regimes,  were  2.  4  x  ]01Z  and  2  x  loll  respectively.  The 
effect  of  cumulation  in  the  cited  experiments  is  concluded  to  be  a  minor  one. 


3S 


Target 
radius,  cm 

Laser 
energy,  j 

Mean 

[Neutron  output  per  pulse 

temp. ,  ev 

exp. 

calc. 

2,50. 10- 2 

600 

40 

- 

- 

1,25«  10-2 

202 

120 

• 

102 

5,50*  lO"3 

214 

840 

3*  106 

8*i07 

3,00*  10- J 

232 

4*10J 

- 

!•  10 ,0 

Table  I 


Novikov,  N.  P.  ,  and  A.  A.  Kholodilov.  Destruction 
of  thermoplastics  by  the  combined  action  of  gas  and 
powerful  thermal  flux.  I-FZh,  v.  22,  no,  4,  1972, 
618-626. 


This  paper  is  a  repeated  treatment  of  an  experiment  reported 
by  the.  authors  previously  (Effects  of  High  Power  Lasers,  Dec.  1971,  48), 
in  which  the  destruction  characteristics  of  several  polymers  are  compared 
under  combined  laser  and  hot  gas  impact.  In  the  present  case  only  PMMA 
and  polystyrene  specimens  were  used;  surface  heating  was  provided  by 
a  COz  laser  plus  a  coaxial  high-speed  flow  of  heated  nitrogen  over  the 
specimen.  The  gas  jet  diameter  was  more  than  double  the  specimen  diameter 
so  that  the  process  could  be  treated  as  one-dimensional  and  stationary.  The 
resulting  liquefaction,  cavity  formation  and  ejection  rate  of  material  are 
discussed  as  functions  of  beam  power  density  and  jet  velocity;  the  conclusions 
are  as  stated  in  the  cited  earlier  work.  The  main  emphasis  is  on  the  differ¬ 
ences  in  destruction  characteristics  which  depend  on  the  chemical  structure 
of  the  target  material.  Thus  polystyrene  shows  a  monotonic  rise  in  destruc¬ 
tion  rate  with  beam  intensity  and  gas  velocity,  whereas  PMMA  may  show  a 
definite  peak  in  destruction  rate  for  the  same  heating,  as  seen  in  Fig.  1. 

This  is  evidently  caused  by  a  temporary  shielding  effect  by  ejecta  in  PMMA 
for  two  of  the  four  curves  in  Fig.  1(a),  which  dissipates  at  higher  gas  veloci¬ 
ties;  no  similar  effect  was  found  for  polystyrene.  Rough  calculations  were 
also  made  of  the  amount  of  ablated  rraterial  for  the  polystyrene  target,  as 
a  function  of  beam  intensity  and  flow  rate.  An  extensive  theoretical  analysis 
of  the  observed  destruction  mechanisms  is  included. 
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Fig.  1.  Linear  destruction  rate  UL  as  a  function  of 
laser  intensity  and  gas  velocity, 
a-  PMMA 
b-  polystyrene 

Surface  temperatures  are  shown  on  each  curve. 


Krasyuk,  I.  K.,  and  P,  P.  Pashinin,  Breakdown 
in  argon  and  nitrogen  from  a  picosecond  laser 
pulse  at  0.3  5  micron  wavelength.  ZhETF  P,  v.  15, 
no.  8,  1572,  471-473. 

Optical  breakdown  triggered  in  Ar  and  N?  by  the  second  harmonic 
emission  from  a  ruby  laser  was  studied  to  ascertain  the  breakdown  mechan¬ 
ism  from  30-50  psec.  pulses  at  0.  35p.  wavelength.  Breakdown  threshold 
*th  was  measured  in  the  gases  at  a  pressure  in  the  400-4500  torr  range  in  an 
experimental  arrangement  analogous  to  one  described  by  the  authors  and 
A.  M.  Prokhorov  (ZhETF  P,  v.  9,  1969,  581).  The  power  of  the  filtered 
second  harmonic  emission  was  measured  with  a  resolution  equal  to  or  better 
than  20  psec.  The  emission  peak  corresponding  to  the  limit  of  visibility  was 
assumed  to  be  1^.  The  experimental  plots  (Fig.  1)  indicate  that  breakdown 
is  triggered  by  multi-photon  ionization  of  gas  atoms  or  molecules.  This 
mechanism  is  confirmed  by  experimental  data  obtained  by  the  authors  and 
A.  M.  Prokhorov  (ZhETF,  v.  58,  1970,  1 606)  at  0.  69P-  wavelength.  Analysis 
of  the  cited  data  and  that  of  other  authors  reveals  that  the  quasiclassic  formula 
derived  by  Keldysh  adequately  describes  the  relative  decrease  in  IK  ,  i.e.  , 
the  increased  probability  of  photoionization,  with  the  increases  in  frequency 
ot  optical  emission.  In  contrast,  no  theory  exists  to  explain  the  fact  that 
It h  in  ^r  a,,,l  al-so  decreases  when  the  breakdown  is  triggered  by  a  20  nsec 
laser  pulse  at  the  0.  35|j.  wavelength;  accepted  avalanching  theory  predicts 
instead  a  monotonic  rise  in  Ij.^  with  laser  frequency  in  the  nanosecond  case. 
Hence  different  breakdown  mechanisms  must  be  considered  in  the  picosecond 
and  nanosecond  pulse  cases. 


I 
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Fig.  1.  Experimental  plots  of  1^  vs.  pressure: 
o-  in  argon,  A  -  in  nitrogen 


Kostylev,  V.  M. ,  and  N.  V.  Komarovskaya.  Energy 
transfer  in  a  medium  of  low  optical  density.  I-FZh, 
v.  'Ll,  no.  5,  1972,  907-912. 


An  experimental  study  was  made  of  the  radiative  energy  transfer 
in  optically  thin  loose  fibrous  layers  bound  by  diffusely  radiant  and  reflecting 
surfaces.  Allowance  was  made  for  the  effects  of  induced  radiation  and  scat¬ 
tering  from  the  medium  in  approximation  of  local  thermodynamic  equilibrium. 

The  fibrous  layers  were  made  of  a  superthin  (l-2p)  fiberglass  or~30p  thick 
caprone  fibers  bound  by  oxidized  aluminum  and  copper  or  polished  aluminum 
surfaces.  The  effective  thermal  conductivity  X  T  of  the  plane-parallel 
optically  thin  layers  was  measured  in  an  electric  calorimeter  with  a  special 
heat-insulating  shield  in  high  vacuum.  The  maximum  XT  error  was  5%. 

The  experimental  XT  data  are  plotted  in  Fig.  1  in  comparison  with  the  theoret¬ 
ical  Xt(t)  dependence  calculated  from 


X  At^T’L 


(1) 


where  X  is  the  radiative  thermal  conductivity  of  an  optically  dense  layer,  a 
is  the  Stefan- Boltzmann  constant,  T  is  the  arithmetic  mean  of  the  layer 
temperature,  L,  is  the  geometric  thickness  of  the  layer,  and  <r  is  the  reduced 
emissivity  of  the  boundary  surfaces,  which  was  experimentally  measured 
in  the  absence  of  the  loose  fibrous  layer.  Allowance  was  made,  when  calcu¬ 
lating  XT,  for  the  coefficient  71  =  3/2  of  angular  distribution  of  radiation  flux 
intensity  incident  on  the  boundaries.  The  experimental  and  theoretical  XT/X 
versus  r  plots  for  two  different  *  -  values  were  similar.  The  data  indicates 
that  X  dependence  on  L  and  r  of  the  layers  with  diffusely  reflecting  surfaces 
is  de Sv:ribed  with  a  good  approximation  by  (1)  and  similar  formulas.  In 
contrast,  the  experiments  with  a  polished  aluminu  *  boundary  (cold)  surface 
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Fig.  1.  \T  in  W/m/degree  vs.  optical  thickness  r  and 
L(m  m)  of  a  loose  fibrous  layer  of  30  kg/m-* 
volume  density,  «  =0.27.  Solid  lines  are 

calculated  by  (1) 

1-  caprone  fiber,  2-  super-thin  glass  fiber 


revealed  a  significant  dis£repancy  with  the  theoretical  data  calculated  bv  (1) 
It  was  concluded  that  the  p  =  3/2  value  is  acceptable  only  for  gray,  diffusing’ 
boundary  surfaces  of  thin  layers  of  any  r,  i.e.  the  boundary' boundary  has 
the  effect  of  increasing  r  by  a  constant  value. 


Burakov,  V.  S.  ,  P.  A.  Naumenkov,  V.  P.  Ivanov, 
and  G.  A.  Kolosovskiy.  Study  of  the  passage  of 
powerful  laser  radiation  through  an  optically  dense 
plasma.  ZhPS,  v.  16,  no.  2,  1972,  239-242. 

Some  nonlinear  absorption  characteristics  of  laser  propagation 
through  a  plasma  are  described.  The  plasma  used  was  optically  dense 
•  '7/uCm)  anc^  at  ~^ev  *n  a  textolite  capillary  2.9  mm  in  diameter.  Trans- 
missibihty  was  measured  with  a  passively  Q-switched  ruby  laser  generating 
a  30  ns  pulse  at  10  --10  w/cm  ,  A  pronounced  bleaching  peak  was  found 
at  about  10  w/cm  .  as  seen  in  Fig.  1.  The  left-most  point  of  the  extrapolated 
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portion  was  obtained  with  a  free- running  ruby;  it  was  not  possible  to  excite 
he  plasma  in  the  Q-switched  mode  below  10°  w/cm2.  The  absorption  char- 

n«aernStlC  J.8’  /ei^Perature  given  for  the  HI,  Cl,  OI  and  CII  components. 
Results  indicate  the  ion  and  electron  temperatures  vary  almost  in  synchron¬ 
ism.  A  general  conclusion  is  that  in  a  multicomponent,  highly  ionized 
plasma  of  the  type  tested,  deviation  from  equilibrium  concentration  of  elec¬ 
trons  can  be  caused  by  individual  hard-ionizing  elements;  in  the  present 

®®  thls  was  d“e  to  the  CI1  comnonent.  It  follow.s  that  care  must  be  taken 
to  allow  for  nonlinear  absorption  when  using  a  powerful  laser  for  certain 
plasma  diagnostics. 


Askar'yan,  G.  A.,  E.  Ya.  Gol'ts,  and  T.  G. 
Rakhmanina.  Alteration  of  the  propagation 
and  reflection  of  ultrasound  under  the  effect  of 
an  intense  light  on  the  surface  of  a  body  in 

ZhETF ,  v.  62,  no.  3,  1972,  1072-1074. 


A  new  effect  was  investigated  experimentally  in  which  the  propa¬ 
gation  of  sound  is  changed  due  to  intense  light  acting  on  the  medium.  A  flash  of 
an  unfocused  and  unmodulated  neodymium  laser  beam  sharply  reduced  the  reflec¬ 
tion  and  transmission  of  ultrasound  through  the  surface  of  a  steel  plate 
immersed  in  water.  Strong  change  is  found  when  the  surface  temperature  T 
is  high  enough  to  form  vapor  or  gas.  An  expression  was  derived  for  determining 
this  temperature.  In  the  experimental  a  rangement  (Fig.  1),  a  laser  beam 


Fig.  1.  Experimental  sketch  for  investigating 
alterations  of  (a)  propagation  and  (b)  reflection 
of  sound  due  to  an  unfocused  laser  beam  on  the 
surface. 

from  1  falls  on  the  surface  of  a  steel  plate  3,  immersed  in  water  in  vessel  6. 

n  ultrasonic  transmitter  2,  a  piezoelectric  element  of  1  cm  radius,  transmits 
a  directional  ultrasonic  wave  at  a  frequency  of  2  MHz  from  generator  5,  in 
such  a  way  that  the  sound  wave  passes  through  the  surface  region,  illuminated 
by  the  laser  pulse.  The  piezoelectric  receiver  2'  records  the  ultrasound 
radiation,  passing  through  (Fig.  la)  or  reflected  by  the  plate  (Fig.  lb).  Two 
series  of  experiments  were  conducted;  a)  using  a  free-running  lasermth  a 
maximum  energy  ~10  joules,  pulse  duration  .  5  msec,  and  beam  radius 
1  cm;  and  b)  using  an  unfocused,  Q-switched  laser  with  a  pulse  width  of 
30  -40  nsec.  Seven  oscillographs  obtained  during  the  experiments  are  given 
to  show  the  effects  of  light  on  the  reflection  of  sound.  The  build-up  time  of 
the  interaction  effect,  connected  with  the  formation  of  vapor-gas  bubbles  or 
of  a  non-uniform  film,  was  commensurate  with  the  energy  release  time. 
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da"a‘ty  -  a,  a  Ugh, 

the  laser  pulse  duration  (  ~1  msec)  At  vr  uSfC'  whlch  Slgnificantly  exceeds8 
a  rapid  film  formation  was  observed  which  of*' **  pov?TS  with  Q'-odulaHon 
the  sound  reflection.  Possible  applications  .eltm“iated  Propagation  and  changed 
elimination  of  reflection  and  prop^gatfon  of  sn,  ^ ^  'he  inter™Pti°n  and  8 
of  sound,  when  short  laser  ‘^ulses^are'usel.801111'’  “d  UltraraP^  modulation 

Zakharov,  V.  P.  ,  y.  N.  Chugayev  V  I  7al- 
Yu.  G.  Poltavfso,,  ci-  j  ,,,  y  ’  v*  Zaliva,  and 
of  thin  carbon  film’s  rTTFr^r—n1  g*~aphitizatirm 
SsEST  Uh'Zh,  n”  ^^l^gSSirCTSiKr- 

dcscribt-d,  which  complemtnL'^hT work^T/akh23''0"  ° ’  8  carbon  film  is 
March  1972  monthly  report  p  9  and fnr / ^kharov  reported  previously 

User  a  type  IFP  flaLlamp  was^hosenlr.h  ”POrl’  P'  ,30>-  Instead  of  a 

peak  a,  0  4  micron,  and  used  to  irradiate  a  lo-’?8'””  Case’  navi"8  a  spectral 
various  distances  and  flash  intensities  u  1-  Cm  PUre  carbon  film  at 

to  increasing  change  in  optical  transmiss  ^  bf  ^  °f.the  exPeriment  was  devoted 

induced  graphitization  process.  A  typical  result  l™’  an  indeX  of  the 

«  typical  result  is  shown  in  Fig.  1. 
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«rton  fifm.an6<!  ”  tranami‘tad  flux  of  irradiated 

a-  after  first  pulse; 
b-  after  second  pulse 
major  divisions  lO-4Sec. 

oM00h0oc  peerap„fi!”  “t,™  }ajTn  ™ ba  raiaad  «°  the  order 
where  graphitization  occurs-  the  red  ,,gLlre  shows  the  inflection  interval 
for  subsequent  pulses.  £ thi, a  g"ohiteSPre  °f  CUrve  (b»  tb“  ba«  true 
Sec.;  correspondingly  increased  tin?  ph  *  tl0n  occurred  in  about  10-5 
treating  rate,”  The'  authors  .u.  J«  a  t  ""'i  re<,Uired  for  reduced 
with  a  rapid  formation  oferystamtes  IndfoUnw1’2!’';,'585  occurri"8-  beginning 
of  crystallite  reorientation  and  groupins  tato  ,h.  m°rC  tended  period 

tesu  „n  r-f  transmissibUity  of  the  eV.ed  film  Sr  ^d ^ 
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Mirkin,  L.  I.  Dynamic  deformation  of  low-carbon 
steel  from  the  effect  of  a  laser  beam.  IN:  Sbornik. 
Vysokoskorostnaya  deformatsiya.  Moskva,  Izd-vo 
Nauka,  1971,  109-112.  (RZhMekh,  8/72,  #  3V1469) 

Structural  effects  are  studied  in  low-carbon  steel  exposed  in 
vacuums  to  focused  laser  pulses  in  the  1  millisecond  range,  at  energies 
up  to  35  j.  The  amount  and  distribution  of  resultant  twinning  was  measured. 
A  physical  model  of  the  beam  action  is  postulated  to  explain  the  simultaneous 
pre  sence  of  a  thermal  and  mechanical  interaction  zone. 


Gurevich,  V.  I.  Pulse  forms  of  a  periodic  point 
source  of  heat  on  the  surface  of  a  large  body. 

FiKhOM,  no.  2,  1972,  19-22. 

A  study  on  the  effect  of  pulse  shape  on  laser  interaction  with 
metals  was  mentioned  by  Baranov,  Gurevich,  and  Heinrichs  in  a  previous 
report  (April  Monthly  Report,  p.  4).  In  the  present  paper  Gurevich  gives 
a  more  extended  analysis  of  pulse  shape  effect.  The  model  assumes  a 
periodic  pulse  from  either  a  stationary  or  moving  source,  and  is  used  to 
calculate  a  limiting  temperature  field  in  the  impact  region  at  the  conclusion 
of  the  laser  pulse;  for  convenience  a  dimensionless  temperature  0  .  is  intro¬ 
duced.  Analytical  expressions  for  0j!  (fixed)  and  0,"  (moving  source)  are 
then  obtained  in  terms  of  beam  parameters  and  the  Fourier  (Fo)  and  Peclet 
(Pe)  criteria.  A  comparison  of  pulse  shape  effect  on  0  •'  is  seen  in  Fig.  1 
for  the  fixed  source  case,  showing  the  maximum  effect  of  a  sawtooth  pulse 


Fig.  1.  Temperature  0  vs.  duty  factor  S 

1-  sawtooth  pulse,  vertical  leading  edge; 

2-  rectangular  pulse; 

3-  sawtooth,  vertical  trailing  edge 
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with  vertical  leading  edge,  at  the  higher  Fo  levels.  It  follows  that  this 
form  would  be  preferable  for  fast  local  heating  of  a  limited  surface  area, 
despite  the  fact  that  a  sin^  pulse  generally  gives  most  efficient  energy 
transfer  in  a  given  pulse  width.  Tabula’ ed  results  are  also  included  com¬ 
paring  the  limit  temperatures  for  a  moving  beam  in  terms  of  the  cited  pulse 
shapes,  again  showing  the  superiority  of  the  vertical-rise  sawtooth. 


Kapel'yan,  S.  N.  ,  and  A.  M.  Yudovin.  Duration 
of  vaporization  after  termination  of  a  powerful 
thermal  flux.  DAN  BSSR,  no.  3,  1972,  214-216. 

Theoretical  expressions  are  developed  which  define  post-pulse 
vaporization  duration,  as  well  as  depth  of  the  vaporization  layer,  for  the 
case  of  laser  irradiated  metals.  The  work  is  based  on  heating  concepts 
reported  by  Anisimov  (Effects  of  High  Power  Lasers,  Dec.  1971,  p.  24) 
and  uses  his  thermophysicai  model.  This  asserts  that  the  thermal  field 
at  the  conclusion  of  a  rectangular  pulse  can  be  given  by 

T  (x)  To  e.\|>  ( —  (hr  i  (1) 

where  1  //3  =  a/vQ  is  a  characteristic  dimension  of  the  heated  region,  and 
T0*  is  the  temperature  at  the  vaporization  front.  Following  the  laser  pulse 
the  vaporization  front  will  continue  to  expand  until  T*  drops  to  vaporization 
threshold.  This  interval  can  be  found  from  a  transcendental  equation 
expressing  vapor  kinematics  and  target  thermal  parameters;  the  authors 
obtained  solutions  by  simple  iteration  using  a  Minsk-22  computer.  Results 
for  several  metals  are  seen  in  Fig.  1,  showing  that  at  a  given  laser  flux 


Fig.  1. 
vs. 


Post -pulse  vaporization  duration 
laser  intensity 
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the  longest  post-vaporization  will  occur  in  Cu  and  Al,  and  the  shortest 
in  Pb.  The  sharp  decrease  in  post-pulse  vaporization  with  rise  in  pulse 
intensity  is  shown  to  be  consistent  with  the  thermal  storage  mechanisms 
in  the  impact  region. 


Pustovalov,  V.  K.  Self-similar  gas  motion 
behind  a  shock  wave  front  sustained  by 
radiation’  DAN  BSSR ,  no.  12,  1971/1079- 
1081. 

The  author  analyzes  a  simple  model  which  relates  to  the  optical 
plasmatron  described  by  Rayzer  and  others  (see  for  example  the  February 
1972  Monthly  Report,  p.  11),  i.  e.  a  continuous  local  plasma  sustained  by 
a  laser  beam.  The  present  model  assumes  a  half-space  x>  0,  filled  with  a 
cold  ideal  gas  of  constant  density  p  ,  where  the  surface  x  =  0  is  the  boun¬ 
dary  between  the  gas  and  vacuum.  At  time  t  —  0  a  strong  shock  wave  begins 
to  propagate  from  the  boundary  in  to  the  gas,  impelling  the  gas  to  expand 
into  the  vacuum.  Energy  from  an  optical  flux  q  is  absorbed  by  the  gas  and 
uniquely  determines  the  propagation  of  the  shock  wave;  gas  expansion  is 
assumed  to  be  adiabatic.  Using  this  model,  the  author  develops  self-similar 
equations  in  Euler  coordinates  defining  gas  pressure,  density,  velocity 
in  terms  of  adiabatic  index  7  and  the  self-similar  index  £  .  It  is  shown 
that  for  the  sustained  shock  condition  7  must  lie  between  1  and  1.  5.  The 
case  of  q<  0,  i.  e.  energy  radiating  from  the  shock  wavefront,  is  also 
briefly  considered. 
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VIII.  PARTICLE  BEAMS 


Gaponov,  V.  A.  ,  and  V.  S.  Nikolayev.  Accelerator 
tube.  Author's  certificate  USSR  no.  299989,  published 
March  2b,  1971,  2  p. 

An  accelerator  tube  with  magnetic  focusing  lenses  is  introduced  . 
The  tube  consists  of  sectionalized  insulated  segments  and  electrodes.  With  the 
aim  of  increasing  the  electric  field  gradient  along  the  tube,  plates  are  placed 
on  the  electrodes  near  the  aperture  for  charged  particle  beam  transmission. 
The  plates  are  located  on  alternate  electrodes  in  a  diametrically  opposite 
manner,  forming  a  periodic  inclined  system..  The  plates  are  comb-shaped 
with  teeth  equivalent  to  elementary  Faraday  cylinders.  Details  of  the 
accelerator  tube  are  shown  in  Fig.  1. 


Fig.  1.  Section  of  accelerator  tube. 

1-  accelerator  ring 

2-  split  metallic  electrodes 

3-  internal  electrodes 

4-  focusing  magnetic  lenses 

5-  apertvires 

6-  flat  plates 
E-  electric  field 
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Kul'man,  V.  G.  ,  E.  A.  Mirochnik,  and  Y.  M. 

Pirozhenko.  Linear  charged  particle  v.ceelerator. 

Author's  certificate  USSR  no.  279822,  published 

March  26,  1971,  2  p. 

A  linear  charged  particle  accelerator  operating  on  a  tr/2  wave 
is  described.  The  accelerator  consists  of  H-profile  accelerating  and  coupling 
elements  with  coupling  windows.  To  increase  the  coefficient  of  cell  coupling, 
the  cells  are  in  the  form  of  ring  resonators.  The  connecting  windows  between 
the  accelerating  and  coupling  cells  arc  placed  with  a  mutual  relative  offset 
in  azimuth.  The  design  provides  for  simple  welding  of  sections  and  ease 
of  evacuation.  Details  are  shown  in  Fig.  1. 


Fig.  1.  Section  of  linear  accelerator. 
1,  2-  Toroidal  accelerating  cells 
3 -  Coupling  cells 
4,  5-  Coupling  windows 


Reproduced  from 
be.<r  available  copy. 
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Rukhadze,  A.  A.,  and  V.  G.  Rukhlin.  Injection  of 
a  relativistic  electron  beam  into  a  plasma.  ZhETF, 
v.  61,  no.  1,  1971,  177-189. 

This  article  deals  with  theoretical  asymptotic  investigations  of 
induced  charges,  currents,  and  electromagnetic  fields,  resulting  from  the 
injection  of  an  electron  beam  into  a  plasma.  The.plasma  was  assumed  to  be 
free  of  other  external  fields.  Only  high-speed  processes  were  taken  into 
account;  thermal  movements  of  particles  were  neglected.  Mathematical 
expressions  were  obtained  for  induced  charges,  currents,  magnetic  fields, 
and  dielectric  permeability.  The  problem  is  represented  by  a  series  of 
mathematical  terms,  so  that  the  general  formulas  obtained  can  be  applied 
to  systems  at  any  moment  of  time  during  the  injection  of  a  finite  length  elec¬ 
tron  beam  into  the  plasma.  At  conditions  rQ>  c/u  ,  where  rQ-  beam  radius 
and  wp-  plasma  frequency,  plasma  perturbations  a*re  localized  in  the  region 
of  the  electron  beam  itself;  and  currents  induced  in  the  plasma  tend  to  com¬ 
pensate  the  beam  magnetic  field,  facilitating  its  injection  into  the  plasma.  If 
the  beam  injection  time  t>  X0  =  y-1  (rQu>  _/c)2  (where  7  is  the  plasma  electron 
collision  frequency),  magnetic  field  compensation  occurs  at  distances 
z^  *"■  zo  =  u r  o  ^rorn  heam  front,  where  u  is  the  directed  electron  velocity. 

I  or  a  high-current  electron  beam  when  the  magnetic  energy'of  beam  current 
exceeds  the  kinetic  energy  of  electrons,  simplified  injection  into  a  dense 
plasma  is  possible  only  for  the  case  where  wp  >  c/r?  and  T<  r Q, 


Vagin,  Yu.  P. ,  G.  L.  Kabanov,  Yu.  A,  Medvedev, 
and  B.  M.  Stepanov.  Method  of  visualizing  space 
distribution  of  dose  in  a  powerful  pulsed  high-speed 
electron  beam.  Atomnaya  energiya,  v.  52,  no. 1, 

1972,  73-75. 

The  scattering  field  of  high-speed  electrons,  and  the  luminosity 
field  in  air  were  investigated  by  photoelectronic  ahd  photographic  methods; 
and  luminosity  field  characteristics  were  compared  with  results  of  approxi¬ 
mation  theory  on  multidimensional  electron  scattering.  Electron  beams  of 
1  and  4  Mev  were  used  with  a  pulsed  electron  current  of  0.  1  and  3,0  a,  a 
pulse  duration  of  2  and  1.2/usec,  and  a  frequency  of  400  and  25  Hz,  respec¬ 
tively.  Electron  beams  were  dispersed  in  air  forming  a  typical  luminous 
cone.  A  photoelectronic  detector  (an  FEU  and  Faraday  cylinder)  simultaneously 
recorded  the  luminous  intensity  and  electron  current.  Luminescence  was 
observed  in  a  small  volume  of  air  (~lcm3),  placed  in  the  Faraday  cylinder  and 
confined  by  an  inlet  diaphragm  (~1  cm).  The  results  are  given  of  measure¬ 
ments  of  the  luminescence  intensity  and  the  electron  current  in  a  lateral 
cross-section  of  an  electron  beam  at  a  distance  of  3  cm  from  the  outlet  window 
of  the  accelerator  (Fig.  1).  Measured  by  three  different  methods, 
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Distance  from  beam  axis,  cm 

Fig.  1.  Lateral  distribution  of  the  intensity  I  of 

luminous  air.  (  *  -  photoelectronic  method; 

O-  photographic  method)  and  electron 
current  with  1  Mev  energy  (A) 

the  results  are  in  good  agreement  and  deviate  from  the  mean  Gaussian  dis¬ 
tribution  curve  by  not  more  than  ±  10%.  The  comparison  of  experimental 
results  with  calculations  using  approximation  theory  on  multidimensional 
electron  scattering  is  plotted  in  Fig.  2.  These  results  agree  well  for  the 
4  Mev  electron  beam;  however,  for  the  1  Mev  beam  satisfactory  agreement 
began  to  appear  only  at  a  distance  of  about  10  cm. 


Fig.  2.  Relationship  of  electron  scattering  angle  0 
in  air  to  distance  L  along  the  beam  axis. 

1,2-  calculated  curves  for  1  and  4  Mev  electrons 
•  ,A-  experimental  curves  for  1  and  4  Mev 
electrons 


52 


IX.  ELECTROMAGNETIC  FIELD  INTERACTION  WITH  MATERIALS 


Balkarey,  Yu.  I.  and  E.  M.  Epshteyn. 
Coulomb  shielding  in  a  strong  electromagnetic 
field.  FTT,  no.  3,  1972,  741-745. 


The  effect  is  analyzed  of  a  strong  electromagnetic  field  on  the 
shielding  of  a  static  charge  by  an  electron  cloud  in  the  electron  plasma  of  a 
crystal.  A  dipole  approximation  is  introduced  on  the  assumption  that  the 
e-m  wave  amplitude  varies  over  a  distance  greater  than  the  charge  dimension, 
the  Debye  shielding  radius,  and  the  electron  oscillatory  amplitude.  In  a 
collisionless  plasma,  the  condition  o)p  t  >  1  must  be  satisfied,  where  io  is 
the  plasma  electron  frequency  and  x  is  the  electron  relaxation  lifetime.  PAn 
integral  equation  of  the  scalar  potential  function  V>(q,t)  was  derived  in  an 
approximation  of  chaotic  phases  from  the  wave  equation  using  Poisson's 
equation.  The  potential  describes  the  static  charge  field.  Using  the  known 

expansion  of  the  Bessel  function  Js  (z),  the  constant  component  of  V’is  expressed 
by 


?«(«•)=  j  fofa)  3 


J  (-"P  ?-  2j  t  (q,  sU) 


(1) 


■ 

where  r  is  the  distance  from  the  shielded  charge  to  the  observation  point, 
q  is  a  coordinate  in  the  Hamiltonian  of  the  analyzed  system,  p  (q)  is  the 
Fourier  component  of  a  static  charge,  a  is  the  oscillatory  amplitude  of  electrons 
in  the  wave  field,  e  (q)  is  the  static  dielectric  constant,  and  to  is  the  electro¬ 
magnetic  field  frequency.  By  introducing 


1 

fa) 


(nq) 
(q.  •«) 


(2) 


in  (1),  a  formula  for  the  effective  static  dielectric  constant  €  eff (q)  is  obtained. 
The  formula  (1)  together  with  (2)  describe  the  effect  of  an  electromagnetic 
wave  on  static  distribution  of  . 


The  component  <PQ  (r)  of  <P  is  calculated  on  the  assumptions  that 
the  shielded  charge  is  a  point  charge  (  p  (q)  =  Ze),  to  >  <op,  and  r  satisfies 
inequalities  r  >  \  and  Xr  Si,  where  \  is  the  characteristic  electron  wavelength 
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and  X  is  the  Debye  radius  reciprocal.  Under  these  assumptions, 


4r.Ze  | 

•  dq  c‘*'r 

x2  1 

f«>  #  ,  I 

CQ  * 

I  (2r.)J  qt  L1 

ql  -j-  *2  K  (a*l)J 

(1)  becomes 
(3) 


When  a  is  much  smaller  than  the  Debye  radius  (  X  a  a  F  <  1),  the  formula 
(3)  is  approximated  by 


.  l  Zci  ( c~“  f  F  2  /  /  1  IN  s-1 

Vo  (r)  = —  \~7T  L1  +  ~  (cos- 3  -1-  ( T  +  W)  <3  cos2  a  -  I  ))J  - 

F'l  1 

—  ^i(3c°S2a-  nj  , 

where  R  =  X  r,  and  a  is  the  angle  between  the  r  and  EQ  vectors. 


(4) 


In  the  presence  of  an  electromagnetic  field,  it  follows  from  (4) 
that  the  potential  distribution  is  anisotropic  and  contains  exponentially  decreas- 
ing  terms  which  predominantly  contribute  to  the  f  Q  value  at  long  distances. 

An  electromagnetic  wave  thus  induces  a  "shielding  breakdown".  At  R  — - >  oo  , 
acquires  a  quadrupole  instead  of  a  dipole  form,  because  the  average 
dipole  moment  is  zero  in  a  high-frequency  field.  If  R  >  max  (2F,  1)  and  a=  0, 
it  follows  from  (4)  that  V>Q  acquires  the  asymptotic  form  <f> Q  =  -(ZeX/t  o) 
(F2/r3),  regardless  of  the  F  value. 


Analogous  mathematical  operations,  using  linearized  equations 
of  motion  and  continuity,  led  to  the  conclusion  that  formulas  (1)  and  (2)  are 
also  applicable  for  frequent  collisions  in  an  electronic  gas,  when  (o  I  £  1. 

F ormulas  (3)  and  (4)  are  also  satisfied  if  the  condition  ;>  <uc  (where  <or  = 
D  X  2  is  the  reciprocal  of  the  Maxwellian  relaxation  lifetime  and  D  is  the 
diffusion  coefficient)  is  substituted  for  r  >  \  and  Xr  2:1.  In  contrast  to  the 
collisionless  plasma,  the  formulas  of  a  and  e  (q,  (o  )  are  explicit. 


In  summary,  "breakdown"  of  shielding  occurs  at  a  sufficiently 
high  e-m  field  frequency  regardless  of  the  status  of  collisions  in  the  electron 
gas.  Numerical  calculations  show  for  example  that  the  condition  fi  ^  a>  p  is 
satisfied  for  n  -  InSb  with  an  average  electron  concentration  N0=  10^  cm-3  at 
fl  =  1014  sec-*.  At  low  temperatures,  EQ  =  105  V/cm  is  required  to 
achieve  F  =  l.  The  condition  fl  ^  o> c  is  satisfied  for  n-Si  with  p  ~  10^  ohm  x  cm 
at  =  1010  sec-1  and  room  temperature;  E0  =  1()3  V/cm  is  sufficient  to 
achieve  F=l.  The  examples  demonstrate  that  the  "breakdown"  of  shielding 
can  occur  at  reasonably  attainable  field  values  in  the  i-r  and  shf  ranges. 
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Bozhkov,  A.  I.  ,  and  F.  V.  Bunkin.  Optical 
excitation  of  surface  waves  in  transparent 
condensed  media.  ZhETF,  v,  61,  no.  6, 

1971,  2279-2286. 

The  excitation  of  surface  waves  in  condensed  media  by  two 
mutually  interfering  plane  monochromatic  waves  is  discussed.  Optical 
excitation  of  surface  waves  in  such  a  case  proceeds  by  a  striction  mechanism 
which  involves  a  ponderomotive  force  jump  from  the  normal  to  the  surface 
component.  The  mechanism  is  not  related  to  radiation  absorption  and  hence 
is  primarily  applicable  to  transparent  media.  The  proposed  method  of  surface 
wave  excitation  is  examined  for  two  coherent  monochromatic  waves  incident 
on  a  liquid  surface.  Mathematical  treatment  of  an  equation  of  motion  of  the 
liquid -atmosphere  interface,  with  allowance  for  ponderomotive  forces, 
produced  general  formulas  of  the  velocity  component  Vz  and  the  function 
£  (x,  y,t)  describing  surface  deviation  from  the  z  =  0  plane.  These  formulas 
led  to  the  conclusion  that  two  plane  electromagnetic  waves  stimulate  vibrations 
of  the  liquid  surface  with  amplitude  lc0l .  a  wave  vector  q  =  ktj  -  anc* 

a  beat  frequency  12=  -  it)  ^  In  these  expressions,  ktj  and  k ^  are  the 

wave  vector  projections  on  the  liquid  surface,  and  coj,  oj  ^  are  the  frequencies 
of  the  two  optical  beams.  A  formula  was  derived  for  |C0  |  ^  as  a  function  of 
q,  12,  which  is  applicable  also  to  a  total  Fresnel  reflection,  when  the  liquid 
occupies  a  space  z  >  C  (x,  y,t).  Analysis  of  frequency  characteristic  |A|  ~2 
variations  as  a  function  of  12  in  the  regions  of  low  and  high  12  revealed  that 
|  A  |  -2  for  viscous  liquids  decreases  continuously  with  increases  in  12,  and  the 
1^1  "2  characteristic  for  low  viscosity  liquids  is  a  resonance  frequency  in 
the  region  of  high  12.  For  12  =  0,  the  I  4  Ql  the  static  surface  wave  is  inde¬ 
pendent  of  viscosity  Y  and  in  the  cases  of  12=  12  Q  (q)  aqd  >  12  Q  /vq^  the  surface 
wave  travels  with  I  £  Ql  dependent  on  v  ,  surface  tension  and  a,  dielectric 
constant  of  the  liquid,  and  independent  of  v  and  a,  respectively.  Conditions 
are  given  for  the  application  of  the  I  C  Q|  formulas  to  a  laser  excitation  source. 
Examples  of  numerical  evaluations  of  the  optical  beam  intensities  required  to 
stimulate  surface  waves  in  a  high  or  a  low  viscosity  liquid  show  that  striction 
can  stimulate  these  waves  with  an  amplitude  much  higher  than  thermal  excitation, 
and  in  certain  cases  may  determine  the  radiation  resistance  of  transparent 
laser  materials. 
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Myshenkov,  V.  I.  and  Yu.  P.  Rayzer. 

Ionization  wave  propagating  as  a  result  of 
resonance  quanta  diffusion  and  sustained  by 
shf  radiation.  ZhETF,  v.  61,  no.  5,  1971, 

1882-1890. 

Sustained  propagation  of  an  ionization  wave  in  noble  gases  in  an 
shf  electric  field  E  was  analyzed.  The  wave  propagates  at  a  power  density 
far  below  the  threshold  density  of  gas  breakdown.  The  simplified  scheme 
of  steady- state  propagation  involves:  energy  transfer  from  the  shf  field  to 
electrons,  atom  excitation  to  a  single  resonance  state,  and  ionization  of  the 
excited  atoms  by  electron  shock  and  diffusion  of  plasma  resonance  radiation 
into  undisturbed  layers.  In  the  one-dimensional  model  of  the  reference 
coordinate  system  (Fig,  1)  the  ionization  wave  is  at  rest  in  the  system.  The 
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Fig.  1.  Schematic  parameters  distribution 
in  a  two-dimensional  steady- state  ionization 
wave:  a  -  density  N*  of  excited  atoms;  b  - 
electron  density  N;  c  -  mean  square  shf  field 
E2.  Arrows  indicate  directions  of  propagation 
velocity  u  of  nonionized  gas  in  the  wave  and 
energy  flux  SQ  of  the  incident  electromagnetic  wave. 


wave  propagation  is  described  by  the  simplified  equations: 


udN  /  dx  =  a.Y.Y*,  a  =  va;  {/■), 


(1) 


where  a  is  the  ionization  rate  constant,  V  is  the  electron  velocity,  and 
0  is  the  ionization  cross-section  of  an  excited  atom; 
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where 

D  is  the  diffusion  coefficient,  T  is  the  average  lifetime 

of  the  excited 

state. 

and  I*  is  the  excitation  potential;  and 

clx 2  c*  \  0)  7 

(3) 

where  <*> is  the  electric  field  frequency,  and  e1  is  the  dielectric  constant. 

Approximate  solution  of  (1),  (2),  and  (3),  with  allowance  for 
boundary  conditions,  gives  the  unknown  functions 

N(x)  =  A'oCv(,),  y  =  —  j  N-clx, 

(4) 

where 
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(5) 
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where 

4 - =-^Tfyi+(u*/u)l±ll,  a*  =  iWIT. 

(7) 

and 

~r~=  —  oic1  = 

ax 

(8) 

where 

X 

S  =  t(x)  ==  j(p  —  |u)f/x, 

oo 

(9) 

S  and  Si  =  S  (1  -  p)  are  the  total  and  dissipative  electromagnetic  energy- 
fluxes,  P  and/* are  the  coefficients  of  reflection  and  absorption,  x  is  the 
relative  lifetime  of  an  excited  atom,  and  u#  is  the  characteristic  velocity. 
Using  (4),  (5),  and  (7),  formula  (10)  is  derived, 
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)] 


(10) 


where  &  and  b  are  constants.  The  maximum  electron  density  Nk  in  the 
plasma  is  consequently  related  directly  to  u.  Using  (4),  (5),  and  (8)  and 
approximating  S(x)  by  the  step  function  S  =  Si  at  x<  0,  S  =  0  at  x  >  0,  an 
equation  was  also  derived  which  correlates  u  with  Sl.  Nk  (S^  can  then  be 
calculated  from  (10)  and  finally,  u(SD)  and  Nk(S0).  The  wave  propagation 
velocity  is  given  by 
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where  Y c  =  const,  if  u  <■  u*.  The  formula  for  maximum  density  of  excited 
atoms 


is  derived  from  (6). 

The  formulas  (4)  -  (12)  and  experimental  and  theoretical  data 
from  the  literature  were  used  to  calculate  u,  Nfc,  N^,  Sq,  and  Sj  for  Xe 
at  p  =  3  torr,  under  the  experimental  conditions  described  by  Bethke  and 
Ruess  [Phys.  Fluids  12,  1969,  822].  The  tabulated  data  show  that  the 
calculated  u  values  increase  with  an  increase  in  shf  power,  in  agreement 
with  the  experiment,  but  are  4--7  times  lower  than  the  experimental  u. 

This  discrepancy  is  possibly  connected  with  the  low  a  value  used  in  calculations. 
The  existence  threshold  of  the  ionization  wave  S  =  cE2/4w  was  calculated 
as  0.4  w/cm^  for  Xe  and  1.2  w/cm^  for  Ar,  which  is  in  reasonable  agreement 
with  the  experiment. 


